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,  ABSTRACT 

The  near  wall  mechanics  of  three-dimensional  turbulent  boundary  layers  were 
studied  using  direct  numerical  simulation.  A  fully  developed  turbulent  channel  flow 
was  suddenly  subjected  to  a  transverse  pressure  gradient.  An  explanation  for  the 
(  reduction  in  the  Reynolds  stresses,  observed  in  experiments,  was  provided  in  terms 

of  the  dynamics  of  near  wall  streamwise  vortices.  A  mechanism  for  the  generation  of 
near-wall  streamwise  vortices  in  turbulent  boundary  layers  was  described.  It  was 
shown  that  the  effect  of  three-dimensionality  is  to  change  the  trajectory  of  fluid 
particles  in  the  vicinity  of  these  vortices,  and  how  this  change  leads  to  reductions 
of  turbulence  stresses.  The  proposed  mechanisms  are  consistent  with  numerous 
I  statistical  data  and  correlations  which  are  presented. 


IV 


ACKNOWLEDGEMENTS 


This  work  was  supported  by  the  United  States  Air  Force  Office  of  Scientific  Research 
unders  Grants  87-0285  and  91-0174,  with  Dr.  James  McMichaei  as  the  technical  moni¬ 
tor.  Computer  time  was  provided  by  NASA-Ames  Research  Center.  We  are  grateful  to 
Professor  Javier  Jimenez  for  numerous  stimulating  discussions,  and  to  Professors  Peter 
Bradshaw  and  John  Eaton  for  useful  discussions  and  helpful  comments  on  a  draft  of  this 
report. 


CONTENTS 


Page 

Abstract . iv 

Acknowledgments .  v 

Contents . vi 

Nomenclature .  viii 


Chapter 


1  Introduction  .  1 

1.1  Motivation .  1 

1.2  Objectives .  2 

1.3  Approach  .  2 

1.4  Three-dimensional  turbulent  boundary  layers .  3 

1.4.1  Definition .  3 

1.4.2  Mechanisms  for  generation  of  3D  TBLs  .  3 

1.4.3  Review  of  previous  results .  5 

1.4.4  Structural  “theories’1  for  the  observed  3D  effects . 11 

2  Numerical  and  analytical  considerations . 15 

2.1  Geometry  and  coordinates . 15 

2.2  Numerical  method . 16 

2.3  Parameters  and  grid . 16 

2.4  The  minimal  flow  unit . 17 

2.5  Preliminary  considerations  for  the  3D  turbulent  channel  flow . 18 

2.6  3D  channel  flow  vs  3D  boundary  layers . 20 

2.7  Methods  of  identifying  the  mechanics  underlying  the  3D  effects  ....  21 

3  Near-wall  turbulence  structures  in  2D  boundary  layers . 25 

3.1  Quasi-streamwise  vortices . 25 

3.1.1  Average  size  and  location . 26 

3.1.2  Detection  of  vortices . 26 

3.1.3  Dynamics  of  quasi-streamwise  vortices  in  TBLs  . 27 

3.2  High-  and  low-speed  streaks . 28 

3.3  Near-wall  shear  layers . 30 

4  Results  and  discussions . 43 

4.1  Reynolds  stresses  during  the  initial  period . 44 

4.1.1  Probability  density  functions,  octant-  and  quadrant  analysis  .  45 

4.2.2  The  companion  flows  method . 47 

4.2  Reynolds  stresses  during  the  later  period  of  three-dimensionality  ...  50 

4.2.1  Changes  in  the  Reynolds  stresses . 50 

4.2.2  Identification  of  the  most  affected  events  in  the  late  period  .  51 

4.2.3  Underlying  mechanisms  in  the  late  period  . 53 

4.3  Mean  velocities  . 55 


VI 


Chapter  Page 

4.4  Vorticity  intensities  . 56 

4.5  Two-point  correlations  and  energy  spectra . 57 

4.6  Comparison  with  previously  suggested  mechanisms . 58 

5  Conclusions . 113 

References .  115 


VII 


NOMENCLATURE 


Roman  Symbols 

al 

\J  u'v'2  -f  v'w'2l(^  +  v®  +  ~ufi)  (structure  parameter). 

ft 

2-component  external  forces. 

kx,  kz 

Components  of  wave  vector. 

l 

shift  of  the  structures  in  the  spanwise  direction. 

Lx,  Ly,  Lz 

Dimensions  of  the  comutational  domain  in  the  x,  y  and  z 

direction. 

NX,Ny,Nz 

Number  of  grid-points  in  the  x,  y  and  2  direction. 

p,p 

Pressure  and  averaged  pressure. 

Re 

Reynolds  number.  R e—Ui2JiA_ 

Rei 

Reynolds  number.  Re^  = 

Rij 

Two-point  velocity  correlation  tensor. 

t 

Time. 

t+ 

W«r  2D- 

uT 

Wall  shear  velocity  \JV^I P- 

u(x,z) 

u(x,2  -f  l). 

u,v,w. 

Ur  2D 

Wall  shear  velocity  in  2D. 

u,v,w 

Velocity  components  in  the  x-,  y-  and  2-direction. 

Ui,U2,U3 

U,  V,  W. 

Uoo 

Free-stream  mean  velocity  in  x-direction. 

UL 

Centerline  velocity  of  the  laminar  flow  with  the  same  mass 

flux. 

U,V,W 

Mean  velocity  components  in  the  x-,  y-  and  2-direction. 

Wy,Wy 

dw  j  dW 

W  and  ~5y- 

Woo 

w  in  the  free-stream. 

w 

w  -  Woo. 

W, 

Wall  motion  in  the  2-direction. 

Woo 

Free-stream  mean  velocity  in  2-direction. 

Axis  in  the  initial  mean  flow  direction. 


x 


y 


yto 


z 

z 


Axis  normal  to  the  wall. 


ury/v. 
uT2Dyl  V- 


Axis  normal  to  x  and 
„  ,  i dpf 2 

2  +  2-^7 1  . 


y- 


Greek  Symbols 

a 
7  i 
7 
7s 
7  r 
v 

V2 

£ 


tan(0fl). 

Angle  of  maximum  streamwise  intensity,  see  eqn.  2. Id. 
Mean  velocity  gradient  angle,  see  eqn.  2.1c. 

Mean  velocity  angle,  see  eqn.  2.1a. 

Reynolds  shear  stress  angle,  see  eqn.  2.1b. 


Kinematic  viscosity. 

d2  ,  a2  ,  a2 

+  ^  4  517- 

Channel  half  width. 


n* 

&R 

6x,Sy,6z 


vorticity  components  in  the  x-,  y-  and  z-direction. 
Mean  uix. 

Angle  of  maximum  auto-correlation. 

Grid  dimensions. 


Special  Symbols 

0* 

0 

O' 


quantity  in  rotated  coordinates, 
ensemble- averaged  quantity, 
fluctuating  quantity. 
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CHAPTER  1 


INTRODUCTION 


1.1  Motivation 

Turbulent  boundary  layers  (TBLs)  of  practical  interest  often  have  a  three- 
dimensional  (3D)  mean  velocity  profile.  This  is  usually  the  case  when  the  geometry 
is  complex,  but  there  are  also  many  important  simple  geometries  with  a  3D  TBL. 
Examples  of  the  latter  are  the  flow  over  a  swept  wing  and  the  cooling  of  a  vertical 
wall  by  a  non-vertical  flow  along  the  wall.  The  number  of  detailed  experiments  on  3D 
TBLs  is  still  far  lower  than  in  two  dimensional  (2D)  boundary  layers,  in  part  because 
of  the  difficulties  in  predicting  2D  turbulent  flows,  which  are  presumed  to  be  simpler. 
Despite  the  dearth  of  experimental  results  several  common  features  of  3D  TBLs,  some 
unexpected,  are  reported  in  several  experiments  (see  the  reviews  by  Johnston  1976, 
Cebeci,  1984).  These  features  are  outlined  below: 

First,  the  eddy-viscosity  differs  for  the  streamwise  and  the  spanwise  component 

of  the  Reynolds  shear  stress,  which  undermines  the  use  of  models  based  on  a  scalar 

eddy  viscosity.  Second,  when  a  2D  TBL  is  turned  into  a  3D  TBL  by  action  of,  say, 

rn  2  _ 2 

a  transverse  pressure  gradient,  the  total  Reynolds  shear  stress,  y  u'v'  4 -  v'w'  ,  is 
lower  than  in  a  2D  TBL  with  a  similar  streamwise  pressure  gradient.  This  is  counter¬ 
intuitive  because  the  addition  of  a  spanwise  mean  velocity  gradient  will  generate 
spanwise  shear  stress,  —  Finally,  ♦herp  is  usually  a  drop  in  the  ratio  of  the 

total  Reynolds  shear  stress  to  twice  the  turbulent  kinetic  energy,  called  the  structure 
parameter  aj.  This  quantity  has  a  near-constant  value  of  about  0.15  in  the  outer  part 
of  many  2D  TBLs  and  has  therefore  been  used  in  turbulence  modelling. 

The  present  turbulence  models  have  difficulty  predicting  3D  TBLs;  features  such 
as  the  reduction  in  the  Reynolds  shear  stress  are  not  predicted.  Therefore,  an  im¬ 
proved  understanding  of  these  phenomena  would  be  very  useful  in  the  development 
of  better  models.  Since  a\  and  the  total  Reynolds  shear  stress  are  invariant  with 
respect  to  a  coordinate  rotation  in  the  plane  of  the  wall,  the  reduction  must  be  due 
to  changes  in  the  structure  of  the  turbulence  rather  than  simple  rotation  of  the  stress 
tensor. 
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In  the  near-wall  region  of  TBLs,  well  known  contributors  to  the  Reynolds  shear 
stress  and  turbulent  kinetic  energy  are  the  quasi- streamwise  vortices,  pumping  high¬ 
speed  fluid  toward  the  wall  and  low-speed  fluid  away  from  the  wall  (Moin,  1987, 
Robinson,  1991).  It  is  therefore  believed  that  this  process  may  have  been  altered  in 
3D  flows,  and  the  mechanism  of  this  alteration  can  be  found  through  a  study  of  the 
turbulence  structures  in  3D  TBLs.  Such  a  study  may  also  contribute  to  increased 
understanding  of  the  role  of  these  structures  in  2D  TBLs  as  well. 

1.2  Objectives 

The  objective  of  this  study  is  to  examine  the  mechanics  underlying  the  ob¬ 
served  3D  phenomena  described  above,  using  direct  numerical  simulation.  Turbu¬ 
lence  statistics  will  be  accumulated  in  order  to  see  how  they  are  affected  by  the  three- 
dimensionality.  Then,  the  effects  of  three-dimensionality  on  the  turbulence  structures 
will  be  examined  to  understand  the  mechanisms  leading  to  the  observed  changes  in 
the  turbulence  statistics. 

The  advantage  of  direct  numerical  simulation  is  that  all  essential  scales  are  re¬ 
solved  so  that  modelling  is  not  necessary.  It  is  possible  to  study  the  unsteady  three- 
dimensional  turbulence  structures  directly,  as  well  as  using  conditional  averaging  and 
other  techniques. 

1.3  Approach 

The  computed  flow  that  is  the  subject  of  this  study  is  a  2D  fully  developed, 
pressure-driven,  turbulent  channel  flow  which  is  turned  by  a  suddenly  imposed  span- 
wise  pressure  gradient.  A  schematic  of  the  channel  with  the  coordinate  system  is 
shown  in  Figure  1.1.  The  x-  and  z-axes  are  in  the  directions  of  the  initial  mean  flow 
and  the  initial  spanwise  direction  respectively,  and  the  y- axis  is  perpendicular  to  the 
v.TiIIs.  In  this  flow  all  effects  due  to  three-dimensionalitv  will  grow  outward  from  the 
wall.  It  therefore  differs  from  spatially  developing  TBLs,  such  as  that  on  a  swept 
wing,  where  three-dimensionality  is  generated  in  the  outer  part  as  well  (see  §1.4.2). 
This  work  is  therefore  a  study  of  the  effects  of  three-dimensionality  on  the  near-wall 
turbulence.  Although  the  presently  simulated  flow  is  a  channel  flow,  the  results  are 
relevant  to  turbulent  boundary  layers  because  the  near-wall  features  of  both  flows  are 
known  to  be  nearly  identical. 
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1.4  Three-dimensional  turbulent  boundary  layers 

1.4.1  Definition 

A  turbulent  bounuary  layer  is  three-dimensional  when,  in  coordinates  fixed  to 
the  wall,  the  directions  of  the  projections  of  the  mean  velocity  vectors  onto  planes 
parallel  to  the  wall  vary  with  the  distance  from  the  wall.  It  follows  that  in  3D  TBLs 
the  \f  ocity  gradient  angle,  7 g  =  tan- /■££■),  var>cs  w>lh  the  distance  from  the 
wall.  U  and  W  are  the  streamwise  and  spanwise  mean  velocity  components  in  planes 
parallel  to  the  wall.  The  maximum  change  in  the  velocity  gradient  angle  across  the 
TBL  can  be  used  as  an  indicator  of  the  degree  of  three-dimensionality  and  is  invariant 
with  respect  to  coordinate  rotation  in  the  plane  parallel  to  the  wall.  Note  that  the 
mean  velocity  profiles  in  the  streamwise  and  the  spanwise  directions  can  have  different 
thickness;  this  is  the  case  for  the  present  channel  flow. 

1.4.2  Mechanisms  for  generation  of  3D  TBLs 

Since  the  mean  velocity  gradient  angle  is  non-zero  in  3D  TBLs,  the  governing 
equation  for  can  be  used  to  identify  the  different  conditions  that  lead  to  three- 
dimensionality.  When  the  boundary  layer  approximations  apply,  Qx  %  and 
the  better  known  transport  equation  for  the  mean  streamwise  vorticity  can  be  used 
instead  (Bradshaw,  1987).  We  have  chosen  not  to  make  this  approximation.  Differ¬ 
entiating  the  Navier  Stokes  equation  for  W  results  in, 

DWy  _  dU_  V2W  1  92 p  dwdu  d2^7  d2^2  df: 

Dt  dx  y  pdydz  dx  dy  dydx  dy 2  dydz  +  3y 

where  Wy  s  x  and  z  are  the  coordinates  in  planes  parallel  to  the  wall,  y  is  normal 
to  the  wall,  and  fz  is  the  2-component  of  the  externa!  forces.  ((/,  V,  W)  and  (u\  v',  w1) 
are  the  mean  and  fluctuating  velocities  in  the  (x ,y,z)  directions  respectively.  The 
normal  derivative  of  Wy  at  the  wall  is  given  by  the  Navier  Stokes  equation  for  IT, 

dWy  l,ldP  DW  , ,  d2W  d2W 

A  term-by-term  study  identifies  the  following  sources  for  three-dimensionality  : 

1)  A  spanwise  pressure  gradient,  will  generate  a  continuous  flux  of  Wy  at  the 
wall,  causing  three-dimensionality  to  grow  outward  from  the  wall.  can  also 
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generate  three-dimensionality  inside  the  flow  directly,  which  will  make  the  entire 
TBL  three-dimensional  much  faster.  This  occurs  when  the  spanwise  pressure  gra- 
dient  varies  in  x,  which  would  generate  The  term  in  cqn.  (Ill) 

will  then  be  non-zero  across  the  TBL.  This  effect  is  often  called  “inviscid  skewing” 
since  it  corresponds  to  a  generation  of  mean  streamwise  vorticity  by  skewing  the 
mean  spanwise  vorticity,  fl-,  through  the  tilting  term  For  steady  flows,  the 

ensuing  spanwise  mean  velocity  gradient  angle  is  often  well  approximated  with 
the  Squire-Winter- Hawthorne  equation,  =  '^(77)’  where  U  is  the 

mean  velocity  in  the  streamwise  direction  and  W  is  the  rorresponding  spanwise 
component.  This  is  a  linearized  approximation  for  low  turning  angles,  but  exper¬ 
iments  have  shown  that  it  is  often  valid  for  larger  angles  as  well.  The  resulting 
Johnston-  or  hodograph-plot,  where  W  is  plotted  against  [J  in  coordinates  aligned 
with  the  local  free-stream  direction,  will  then  be  almost  triangular  in  shape.  The 
term  is  negligible  for  thin  TBLs,  but  may  play  a  role  when  the  TBL  is  thick. 

2)  If  the  flow  encounters  a  wall  moving  in  a  direction  different  from  the  local  mean 
velocity  gradient  direction,  three-dimensionality  is  generated  at  the  wall  and  will 
grow  outward.  These  flows  are  called  shear-driven  3D  TBLs.  In  the  case  of  a 
constant  wall  velocity,  the  only  sources  of  Wy  are  the  horizontal  viscous  forces 
between  the  moving  and  the  non-moving  part  of  the  wall  (last  two  terms  in  eqn. 
(1.4.2)).  If  the  spanwise  wall  velocity  varies  in  time,  a  continuous  flux  of  Wy  is 
generated  at  the  moving  wjl  through  the  term  Note  that  this  term  generates 
a  flux  of  Wy  at  the  wall  in  the  same  manner  as  a  spanwise  pressure  gradient.  A 
channel  flow  turned  by  a  constant  spanwise  pressure  gradient  is  equivalent  to  a 
channel  flow  turned  by  a  wall  accelerating  in  the  spanwise  direction  at  a  constant 
rate  (see  §2.5). 

3)  Body  forces  can  generate  3D  TBLs.  Gravity  is  usually  not  important  unless  there 
is  a  large  temperature  difference  across  the  boundary  layer  which  would  generate  a 
buoyancy  force.  This  force  can  turn  a  TBL  on  a  non-horizontal  wall  by  generating 
Wy  across  the  TBL  (Siebers,  1983).  The  Coriolis  force  due  to  Earth’s  rotation 
causes  Earth’s  TBL  to  be  three-dimensional  (Coleman  et  at.,  1990). 

4)  For  certain  geometries,  like  non-circular  ducts,  the  Reynolds  stress  terms  in  the 
transport  equation  for  Wy  generate  weak  three-dimensionality. 

5)  Wy  can  enter  through  the  outer  edge,  notably  when  a  boundary  layer  interacts 
with  a  longitudinal  vortex  (Bradshaw,  1987). 
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These  five  ways  in  which  three-dimensionality  is  imposed  on  a  boundary  layer  may 
affect  different  regions  of  the  boundary  layer  and  with  different  intensities.  They  may 
therefore  lead  to  different  3D  effects  because  the  turbulence  structures  are  different 
near  the  wall  and  in  the  outer  layer.  Additional  effects  may  appear  when  the  trnning 
is  strong.  Another  important  distinction  is  between  equilibrium  and  non-equilibrium 
3D  turbulent  boundary  layers.  In  non-equilibrium  3D  TBLs,  a  flow  which  is  initially 
2D  or  nearly  so,  is  turned  by  a  transverse  force  and  the  observed  3D-phenomena 
occur  while  the  turbulence  adjusts  to  the  new  flow  direction.  When  the  3D  TBL  is 
in  equilibrium,  the  turbulence  has  had  sufficient  time  to  adjust.  An  example  of  the 
latter  is  the  flow  on  an  infinitely  large  rotating  disk  (Liitell  k  Eaton,  1992). 

1.4.3  Review  of  previous  results 

Results  from  pressure-driven  and  shear-driven  3D  TBLs  will  be  discussed  below. 
Three-dimensionality  generated  in  other  ways  is  not  included. 

Pressure  —  driven  flows 

Pressure-driven  flows  are  of  great  practical  interest.  Important  examples  are  the 
flow  over  a  swept  wing  and  flow  around  buildings  and  other  obstacles. 

In  the  flow  over  a  swept  wing,  the  incoming  flow  is  turned  by  a  pressure  gradient 
which  is  roughly  perpendicular  to  ttie  leading  edge  in  plan  view.  Except  at  the  very 
front  of  the  wing,  the  pressure  gradient  rs  adverse  and  will  bend  the  flow  away  from 
the  body  of  the  aeroplane.  In  order  to  reduce  the  number  of  degrees  of  freedom, 
swept  wing  experiments  are  set  up  to  approximate  an  infinitely  long  constant-cord 
wing,  such  that  the  flow  is  homogeneous  in  the  direction  of  the  leading  edge. 

Reynolds  stress  data  for  a  3D  TBL  were  first  reported  by  Bradshaw  k  Terrell 
(1969).  This  was  a  flow  recovering  on  a  flat  plate  behind  a  45°  swept  wing.  The 
three-dimensionality  was  weak  with  a  maximum  angle  of  9°  between  the  tunnel  axis 
and  the  surface  streamline.  The  only  observed  3D  effect  was  a  lag  between  the  stress 
angle  tan_1(p*»)  and  the  mean  velocity  gradient  angle  tan~l(^~/|j~)-  A  practical 
consequence  of  this  is  that  the  eddy-viscosity  is  not  a  scalar.  This  prevents  a  simple 
extension  of  eddy-viscosity  turbulence  models  to  apply  to  3D  flows. 

Experiments  in  a  duct  shaped  to  simulate  a  swept  wing  were  done  at  the  National 
Aerospace  Laboratory  (NLR)  in  the  Netherlands  as  a  response  to  the  need  for  more 
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data  to  test  and  develop  3D  turbulence  models.  The  turbulence  statistics  were  re¬ 
ported  by  van  den  Berg  &  Elsenaar  (1972),  Elsenaar  &:  Boelsma  (1974)  and  van  den 
Berg,  Elsenaar,  Lindhout  &;  Wesseling  (1975).  As  in  the  experiment  of  Bradshaw  and 
Terrell,  the  stress  angle  lagged  behind  the  mean  velocity  gradient  angle.  A  new  find¬ 
ing  was  that  the  structure  parameter  ai  =  yfu  'v'^  4-  v'w'^ f(uf2  -f  t/2  +  u/2),  dropped 
below  the  commonly  reported  value  of  0.15  in  2D  flows.  The  drop  in  ai  in  the  outer 
part  of  the  TBL  occurred  before  the  changes  near  the  wall  had  sufficient  time  to  dif¬ 
fuse  outward.  This  indicated  that  the  inviscid  skewing  mechanism  described  earlier 
played  a  role  in  these  flows. 

In  a  Euromech  meeting  in  Trondheim,  Norway,  (Fannelpp  &  Krogstad,  1975)  the 
existing  turbulence  models  were  applied  to  3D  TBLs  for  which  experimental  data 
existed.  During  this  meeting,  known  as  the  “Trondheim-trials”,  it  became  clear  that 
the  turbulence  models  did  not  perform  well  for  3D  flows  like  the  NLR  “wing”. 

Bradshaw  &  Pontikos  (1985)  performed  an  experiment  with  a  similar  geometry  to 
the  NLR  experiments  except  that  the  adverse  pressure  gradient  was  milder  to  avoid 
separation.  They  were  able  to  confirm  the  3D  effects  found  at  NLR.  The  stress  angle 
lagged  behind  the  mean  velocity  gradient  angle,  though  never  more  than  8°,  and  the 
structure  parameter  dropped  below  0.15  across  the  TBL. 

When  the  flow  is  turned  by  an  obstacle,  the  adverse  pressure  gradient  is  often 
stronger  and  the  turning  more  rapid  than  for  swept  wings.  This  will  usually  make 
the  flow  separate  earlier  and  measurements  must  be  made  over  a  shorter  range  in  the 
streamwise  direction.  The  measured  turbulence  statistics  may  therefore  not  have  had 
sufficient  time  to  be  altered  except  near  the  wall. 

Johnston  (1970)  investigated  the  flow  over  a  45°  swept  rectangular  step  with 
a  very  strong  adverse  pressure  gradient.  The  vertical  pressure  gradient  may  have 
affected  the  three-dimensionality  in  this  flow  through  the  term  in  eqn.  (1.4.1). 
Because  of  the  strong  pressure  gradients,  the  flow  turned  as  much  as  35°  over  a 
streamwise  distance  of  only  five  b<  idary  layer  thicknesses.  As  with  the  earlier 
experiments  of  Terrell  and  Bradshaw,  the  Reynolds  shear  stress  angle  lagged  behind 
the  mean  velocity  gradient  angle.  The  difference  in  angle  was  as  high  as  20°,  much 
larger  than  what  is  found  for  swept  wings. 

Anderson  &  Eaton  (1987,  1989)  investigated  the  TBL  around  a  wedge  that  was 
attached  to  the  wall  and  pointed  in  the  upstream  direction.  Two  wedge  angles  were 
investigated,  90°  and  60°.  As  with  Johnston’s  experiment,  the  3D  effects  were  much 
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larger  than  in  the  swept  wing  experiments  and  were  largest  for  the  largest  wedge 
angle.  Values  of  the  structure  parameter  aj  were  as  low  as  0.03  and  the  ratio  of  the 
spanwise  to  streamwise  eddy  viscosity,  in  coordinates  aligned  with  the  local  mean 
flow,  was  as  low  as  0.2.  The  streamwise  evolution  of  aj  was  dominated  by  the  out¬ 
ward  diffusion  of  near-wall  effects  to  a  larger  extent  than  in  swept-wing  experiments. 
A  comparison  of  existing  3D  data  led  to  the  conclusion  that  the  eddy- viscosity  ra¬ 
tio  decreases  with  increasing  rate  of  boundary  layer  turning,  defined  as  the  ratio  of 
boundary  layer  thickness  to  free-stream  streamline  curvature  radius,  while  being  in¬ 
dependent  of  the  difference  in  mean  velocity  gradient  angle  between  the  wall  and 
the  free-stream.  However,  the  structure  parameter  was  found  to  depend  on  both.  A 
significant  drop  in  ai  and  a  low  eddy-viscosity  ratio  was  also  found  by  Dechow  & 
Felsch  (1977),  who  investigated  the  flow  around  a  wing-shaped  obstacle. 

Moin,  Shih,  Driver  &  Mansour  (1990)  investigated  the  statistical  evolution  of  a 
fully  developed  2D  channel  flow  subjected  to  a  sudden  spanwise  pressure  gradient. 
The  structure  of  this  flow  was  investigated  in  the  present  work.  This  flow  is  different 
from  the  previously  described  flows  in  that  all  three  dimensionality  grows  outward 
from  the  wall  with  no  possibility  of  inviscid  skewing.  The  flow  is  well  suited  for  the 
study  of  near-wall  three-dimensionality  because  the  change  in  turbulence  statistics 
from  their  initial  profiles  to  the  final  steady  state  profiles  is  gradual  enough  for  the 
three-dimensional  effects  on  the  statistics  to  be  easily  observed.  The  channel  flow 
displays  a  drop  in  the  total  Reynolds  shear  stress,  the  turbulent  kinetic  energy  and  the 
structure  parameter  aj.  Since  the  magnitude  of  the  total  Reynolds  shear  stress  and 
the  turbulent  kinetic  energy  (normalized  with  the  initial  shear  velocity  and  channel 
half  width)  increase  with  the  Reynolds  number,  they  will  have  higher  values  in  the 
final  steady  state  than  initially.  The  reductions  must  therefore  be  due  to  the  subtle 
effects  of  three-dimensionality.  From  the  Reynolds  stress  budgets,  they  found  that 
the  decrease  in  —  u'v'  and  u'2  were  due  to  reduced  respective  productions,  while  the 
pressure- velocity  term  — caused  the  reduction  in  t;'2.  The  Reynolds  stress  angle 
lagged  behind  the  mean  velocity  gradient  angle,  with  the  largest  difference  at  the 
wall.  The  mean  velocity,  when  scaled  with  the  total  wall  shear  velocity,  shows  a  log- 
region  with  a  slope  lower  than  the  value  of  2.4  usually  found  for  2D  flows.  Note  that 
this  flow  is  nominally  different  from  a  flow  with  an  impulsively  started  wall  motion 
which  would  be  more  closely  related  to  experimentally  studied  flows  along  a  spinning 
cylinder  (see  below). 
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Shear  —  driven  flows 

Experimental  and  computational  investigations  of  flows  where  wall-motion  gener¬ 
ates  three-dimensionality  will  now  be  discussed.  In  the  absence  of  additional  sources 
for  three-dimensionality  all  3D  effects  will  grow  outward  from  the  wall. 

There  have  been  several  experiments  in  zero  pressure  gradient  flows  along  a  cylin¬ 
der  with  a  rotating  part.  Three-dimensionality  is  generated  by  circumferential  viscous 
forces  at  the  wall  between  the  stationary  and  the  rotating  parts.  If  the  rotating  part  is 
long  enough,  the  flow,  in  coordinates  moving  in  the  spanwise  direction  with  the  speed 
of  cylinder  surface,  Ws,  will  become  2D  with  an  angle  tan-1  {Ws/Uoo)  with  respect 
to  the  cylinder  axis.  The  flow  is  then  said  to  be  collateral.  However,  the  Reynolds 
stresses  and  the  axial  wall  stress,  normalized  with  U0 o,  are  higher  in  the  collateral 
flow  due  to  higher  effective  free-stream  velocity,  \/koo  +  Wf,  an<^  the  radial  pressure 
gradient  (Driver  &  Johnston,  1990).  Therefore,  when  the  flow  enters  or  leaves  the 
rotating  part,  an  adjustment  toward  the  new  2D  profiles  will  grow  outward  as  the 
flow  progresses  downstream.  Results  from  three  cylinder-experiments  are  described 
below. 

Bissonnette  and  Mellor  (1973)  performed  measurements  on  a  cylinder  with  fixed 
forebody  and  rotating  afterbody.  They  found  that  the  flow  remained  collateral  near 
the  wall  when  the  flow  entered  the  rotating  part,  while  further  out  the  flow  was  3D 
and  the  Reynolds  shear  stress  angle  lagged  behind  the  mean  velocity  gradient  angle. 
The  eddy-viscosity  in  the  spanwise  direction  was  in  most  cases  between  0.5  and  0.8 
times  the  streamwise  (axial)  eddy-viscosity.  Using  scaling  arguments,  they  argued 
that  this  was  at  least  in  part  due  to  curvature  effects.  They  also  investigated  the  law- 
of-the-wall  for  this  flow.  Following  a  suggestion  by  Clauser,  they  chose  a  coordinate 
system  that  moved  in  the  spanwise  direction  with  the  wall  and  aligned  with  the  free 
stream  velocity  +  W}.  The  slope  in  the  log-region  was  considerably  lower  than 

in  2D  flows.  Again  scaling  arguments  were  used  and  it  was  found  that  the  curvature 
contributed  to  this  result.  The  spanwise  wall  stress  increased  “instantly”  as  the  flow 
entered  the  rotating  part  and  decayed  downstream,  while  the  streamwise  wall  stress 
increased  slightly  in  the  axial  direction. 

Lohmann  (1976)  investigated  the  same  geometry,  but  with  a  stronger  surface  to 
free-stream  velocity  ratio.  The  spanwise  TBL  did  not  grow  to  the  same  thickness  as 
the  streamwise  TBL  within  the  region  measured,  causing  the  streamwise  boundary 
layer  thickness,  6 ,  to  be  unaffected  by  the  spanwise  wall  movement.  When  the  surface 
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velocity  was  increased,  the  growth  of  the  spanwise  TBL  increased  slightly  and  the 
shape  of  the  spanwise  mean  velocity  profile  outside  y/6  =  0.05  was  altered.  As  in 
the  experiment  of  Bissonette  and  Mellor,  the  TBL  was  collateral  near  the  wall  and 
three-dimensional  further  out,  with  the  stress  angle  lagging  the  mean  velocity  gradient 
angle.  The  Reynolds  stresses  increased  rapidly  near  the  wall  when  the  flow  entered  the 
rotating  part  and  the  increase  diffused  outward  as  the  flow  progressed  downstream. 
Any  weakening  of  the  Reynolds  stresses  due  to  3D  effects  were  therefore  difficult  to 
detect.  However,  for  Ws/f/oo=1.45,  the  axial  wall  stress  decreased  slightly  after  the 
rapid  initial  increase  and  then  continued  to  increase  toward  a  new  steady  state,  and 
for  Ws/Uoo-2.2,  it  remained  constant  for  a  while  after  the  immediate  increase,  and 
then  increased  further  to  a  higher  2D  value.  This  indicates  that  the  3D  effects  tend 
to  reduce  the  axial  wall  stress. 

Driver  &  Johnston  (1990)  reported  on  a  similar  flow.  In  their  experiment,  a  col¬ 
lateral  TBL  on  a  spinning  forebody  became  three-dimensional  as  the  flow  entered 
the  stationary  afterbody.  This  makes  measurements  easier,  but  it  creates  the  addi¬ 
tional  effect  that  the  initial  2D  TBL  is  not  aligned  with  the  axial  direction.  As  the 
flow  enters  the  stationary  cylinder  section,  the  Reynolds  stresses  reduce  to  their  new 
steady  state.  However,  unlike  the  normal  intensities,  the  axial  Reynolds  shear  stress, 
— uV,  fell  below  its  final  steady  state  value  and  then  recovered.  The  reduction  is  due 
to  3D  effects  which  also  caused  the  structure  parameter  a\  to  fall  below  its  steady 
state  value  and  then  recover.  Surprisingly,  a\  was  as  high  as  0.17  on  the  rotating 
part.  The  mean  velocity  gradient  angle  relaxed  toward  zero  more  rapidly  than  the 
Reynolds  shear  stress  angle. 

The  main  purpose  of  Driver  and  Johnston’s  work  was  to  study  the  cylinder  flow  in 
an  adverse  axial  pressure  gradient.  Since  there  was  no  spanwise  pressure  gradient,  a 
mechanism  for  three-dimensionality  did  not  exist  in  the  outer  part  of  the  TBL.  Their 
objective  was  to  investigate  properties  of  the  outward- growing  three-dimensionality 
in  a  TBL  with  an  adverse  pressure  gradient.  Earlier  experiments  on  this  flow  (except 
that  the  flow  progressed  from  the  stationary  part  to  the  rotating  part)  were  performed 
by  Furuya  et  ai  (1977)  and  Driver  &  Hebbar  (1989).  They  found  the  effect  of  the  axial 
pressure  gradient  to  be  minimal.  In  the  absence  of  cylinder  rotation,  the  adverse  axial 
pressure  gradient  tends  to  weaken  the  Reynolds  stresses  near  the  wall  and  increase 
them  further  out.  Driver  and  Johnston  found  that  when  the  flow  is  three-dimensional, 
the  increase  in  the  outer  region  is  less  pronounced.  Along  a  stationary  cylinder  with 
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an  axial  pressure  gradient,  a\  was  weakened  across  the  TBL,  but  its  radial  gradient 
remained  close  to  zero  outside  the  near-wall  region.  In  the  3D  TBL,  aj  did  not  drop 
near  the  wall,  but  further  out  it  dropped  increasingly  with  increasing  distance  from 
the  cylinder. 

Shear  —  driven  flows  in  equilibrium 

Some  studies  have  been  performed  in  the  so-called  equilibrium  flows,  where  the 
turbulence  experiences  the  three-dimensionality  over  a  long  time  and  has  had  time  to 
adjust  to  it.  One  of  the  reasons  for  studying  these  flows  is  that  they  often  have  fewer 
degrees  of  freedom.  The  relevance  of  these  flows  to  the  non-equilibrium  3D  flows  is 
not  yet  clear,  but  they  show  some  non-equilibrium  characteristics.  Three  studies  are 
described  below. 

The  first  of  these  is  the  direct  numerical  simulation  of  a  horizontally  homogeneous 
flow  with  the  free-stream  velocity  vector  rotating  about  a  vertical  axis  at  a  constant 
rate  and  a  fixed  magnitude  (Spalart,  1989).  In  coordinates  rotating  with  the  free- 
stream,  the  flow  is  statistically  steady.  When  the  magnitude  and  rotation  rate  of  the 
free-stream  vector  is  used  for  normalization,  the  Reynolds  number  and  the  distance 
from  the  wall  are  the  only  independent  variables.  The  aj  parameter,  when  plotted 
against  y+,  is  reduced  and  the  amount  of  reduction  increased  with  increasing  Re. 
However,  the  difference  between  the  Reynolds  shear  stress  angle  and  the  mean  velocity 
gradient  angle  was  small  and  changed  sign  several  times  across  the  TBL. 

The  Ekman  layer,  the  prototype  of  planetary  boundary  layers,  was  simulated  by 
Coleman  et  al.  (1990).  In  the  free  stream,  the  pressure  gradient  is  perpendicular  to 
the  flow  direction  and  is  balanced  by  the  Coriolis  force.  The  structure  parameter 
profile  varied  greatly,  but  was  mostly  below  0.15.  It  did  not  seem  to  have  any  region 
with  a  constant  value.  The  difference  between  the  Reynolds  shear  stress  direction 
and  the  mean  velocity  gradient  direction  was  small  in  this  flow  as  well. 

Littell  and  Eaton  (1992)  investigated  a  3D  TBL  on  a  spinning  disk.  Profiles  of 
the  Reynolds  stresses,  normalized  with  the  wall  shear  velocity,  collapsed  well  when 
plotted  against  the  distance  from  the  wall,  normalized  with  the  momentum  thickness, 
for  different  values  of  rotation  rate  and  radius.  Again,  the  lag  between  the  turbulent 
shear  stress  direction  and  the  mean  velocity  direction  was  very  small.  The  structure 
parameter  had  a  surprising  profile  with  a  peak  near  the  wall  and  a  near-linear  drop 
further  out.  Conditional  averaged  two-point  correlations  were  computed  and  based 
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on  these  it  was  suggested  that  a  vortex  generates  weaker  sweep-motion  (t/  >  0,  v'  <  0 
and  -uV  >  2VV^^)  in  a  3D  TBL  if  its  sign  of  rotation  is  opposite  to  the  mean 
streamwise  vorticity  at  the  wall  and  weaker  ejections  («'  <  0,  v'  >0  and  —  uV  > 


1.4.4  Structural  “theories”  for  the  observed  3D  effects 

It  is  known  that  quasi-streamwise  vortical  structures  play  an  important  role  in 
TBLs.  The  weakening  of  the  total  Reynolds  shear  stress  and  other  3D  effects  are 
therefore  believed  to  reflect  the  effect  of  three-dimensionality  on  these  eddies. 

In  the  spinning  cylinder  flow,  Lohmann  (1976)  found  an  increased  contribution 
from  small  scales  and  a  reduced  contribution  from  large  scales  in  the  energy  spectrum. 
He  therefore  suggested  that  the  transverse  shearing  motion  broke  up  the  large  eddies 
into  smaller  eddies  and  that  recovery  represented  the  build-up  of  large  eddies.  The 
increase  in  small  scales  may  in  part  be  an  artifact  of  the  coordinate  system,  since  the 
2D  flow  has  more  small  and  medium  size  scales  in  the  spanwise  direction. 

Bradshaw  &  Pontikos  (1985)  suggested  that  the  eddies  are  tilted  over  in  the 
spanwise  direction  by  the  spanwise  strain.  Since  the  eddies  near  the  wall  have  a  low 
angle  relative  to  the  wall,  while  further  out  they  appear  with  high  inclination,  the 
outer  eddies  are  most  likely  to  be  toppled.  Robinson  (1991)  described  the  dominant 
vortical  structure  in  the  outer  part  of  2D  TBLs,  as  consisting  of  an  arch-shaped  front 
with  high  inclination  angle  and  one  or  two  legs  with  lower  angle  closer  to  the  wall. 
The  strongest  ejections  in  the  outer  part  were  found  inside  the  arch.  These  ejections 
would  be  affected  if  the  arch  is  convected  further  in  the  spanw'ise  direction  than  the 
legs. 

Anderson  &  Eaton  (1989)  suggested  that  near  the  wall  the  mean  streamwise 
vorticity  will  destroy  the  streamwise  vortices  with  opposite  sign  of  rotation  and  roll 
up  into  vortices  with  the  same  sign  of  rotation.  Vortices  with  one  sign  would  dominate, 
which  could  reduce  the  contribution  from  strong  mixing  motion  between  two  vortices 
of  opposite  sign. 

Simpson  &  Devenport  (1990)  suggested  that  fluid  swept  toward  the  wall  will  in 
3D  TBLs  displace  most  low-speed  fluid  in  one  spanwise  direction.  This  fluid,  when 
ejected,  would  generate  positive  u'w'.  However,  they  pointed  out  that  the  fluid  need 
not  eject  outward  to  the  same  level  as  in  a  2D  TBL,  because  the  flow  “has  a  third 
dimension  of  spatial  freedom”.  It  is  not  clear  what  kind  of  mechanism  this  refers  to 
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since  2D  and  3D  flow  terminology  refers  to  the  mean  flow  only  and  all  turbulent  flows 
are  three-dimensional. 

Eaton  (1992)  suggested  that  the  behavior  of  an  imbedded  vortex  in  3D  TBLs 
could  cast  light  on  the  behavior  of  naturally  occurring  near-wall  vortices.  The  latter 
are  much  smaller,  but  are  often  single-legged  with  a  low  inclination  angle.  Shizawa  & 
Eaton  (1990)  introduced  a  single,  large,  longitudinal  vortex  of  either  sign  into  the  3D 
wedge-flow  of  Anderson  &  Eaton  (1989)  and  found  that  vortices  of  both  signs  decayed 
quicker  than  they  would  in  a  2D  TBL.  In  addition,  the  ejection-mechanism  depended 
strongly  on  the  sign  of  the  embedded  vortex.  If  the  vortex  induced  secondary  flow 
near  the  wall  in  the  same  direction  as  the  spanwise  pressure  gradient,  no  secondary 
flow  separation  took  place  and  no  low-speed  streaks  were  generated.  Vortices  of 
opposite  sign  induced  motion  opposing  the  3D  TBL  and  a  strong  outward  motion 
took  place.  Eaton  (1992),  based  on  the  idea  that  ejections  take  place  as  a  result  of 
streaks  becoming  unstable,  suggested  that  the  reduced  number  of  streaks  generated 
less  ejections  and  therefore  less  Reynolds  shear  stress. 

It  should  be  noted  that  these  suggestions  are  not  mutually  exclusive.  The  key 
mechanism  may  depend  on  the  distance  from  the  wall,  the  turning  rate  and  other 
quantities. 

In  Chapter  4  the  results  of  our  study  of  the  three-dimensional  channel  flow  will 
be  described  in  detail.  The  mechanisms  underlying  the  reduced  Reynolds  shear  stress 
will  be  presented  and  compared  with  the  previously  suggested  mechanisms  described 
above. 
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CHAPTER  2 


NUMERICAL  AND  ANALYTICAL  CONSIDERATIONS 

In  this  chapter  a  description  of  the  geometry,  the  numerical  method  and  physical 
and  numerical  parameters  are  presented.  In  addition,  the  tools  used  for  the  analysis 
and  postprocessing  of  the  simulation  data  are  described. 

2.1  Geometry  and  coordinates 

The  simulated  flow  is  initially  a  fully  developed  turbulent  Poiseuille  flow.  This 
is  a  pressure-driven  flow  between  two  infinitely  large,  parallel,  horizontal  walls  with 
a  steady  mean  velocity,  (17(y),0,0).  A  constant  spanwise  pressure  gradient  is  sud¬ 
denly  imposed  which  causes  a  transient,  three  dimensional  evolution  of  the  flow.  The 
perturbed  flow  eventually  becomes  2D  and  statistically  steady  in  a  new  direction.  It 
is  the  transient  behavior  which  is  of  interest  in  this  study.  The  geometry  and  the 
coordinates  for  the  simulated  flow  are  shown  in  Figure  1.1.  The  ar-axis  is  aligned  with 
the  initial  mean  flow  direction,  the  y- axis  is  normal  to  the  walls,  and  the  z-axis  is  in 
the  initial  spanwise  direction.  The  velocity  components  in  the  ( x,y ,  z)  directions  will 
be  denoted  by  (u,v,w),  or  (111,1*2,113)  respectively. 

One  of  the  difficulties  in  studying  three-dimensional  turbulent  boundary  layers 
is  the  choice  of  coordinates.  In  rotated  coordinates  the  velocity  components  will  be 
denoted  as  (u*,  v,  w+).  There  are  several  characteristic  angles,  and  those  used  in 
the  current  work  are  listed  below.  They  will  be  called  the  mean  velocity  angle,  the 
mean  velocity  gradient  angle,  the  Reynolds  shear  stress  angle  and  the  intensity  angle, 
respectively 


7s  =  tan- 1  [TV’/ 1/] 

7  T  —  tan~l[t/u//uV] 

*sdW  ,du, 
7,  =  tan  [g^l 


7 {  —  0.5  tan  l[2u'w'/(ur2  —  u/2)] 


(2.1a) 

(2.16) 

(2.1c) 

(2.  Id) 


The  latter  is  the  angle  of  maximum  streamwise  intensity. 


2.2  Numerical  method 


The  2D  turbulent  channel  flow  was  computed  by  Kim,  Moin  and  Moser  (1987). 
They  used  direct  numerical  simulation  technique  which  resolves  all  the  essential  scales 
and  no  modelling  is  needed.  The  same  computer  code  is  used  here  in  the  simulation 
of  the  3D  channel.  Several  stored  velocity  fields  of  Kim  et  al.  were  used  as  initial 
conditions  for  different  realizations  of  the  3D  flow.  A  brief  description  of  the  numerical 
method  is  given  below.  For  a  detailed  description,  refer  to  Kim  et  al.  (1987). 

The  numerical  method  is  a  spectral  method  with  pressure  eliminated  from  the 
equations.  Fourier  expansions  are  used  for  computation  of  the  horizontal  spatial 
derivatives  and  the  Chebychev  expansion  is  used  for  vertical  derivatives.  The  vis¬ 
cous  terms  are  time-advanced  using  the  Crank-Nicholson  method,  while  the  Adams- 
Bashforth  method  is  used  for  the  non-linear  terms.  The  discretized  equations  are 
Fourier  transformed  in  the  horizontal  plane  and  then  solved  using  the  Chebychev-tau 
method  (Gottlieb  k  Orszag,  1977)  for  each  wave  number.  The  governing  equation 
for  the  vertical  component  of  the  vorticity,  u>y,  is  solved  first,  followed  by  the  fourth 
order  equation  for  V2u  which  is  separated  into  two  second-order  equations.  The  ex¬ 
tra  boundary  conditions  needed  at  the  walls  are  provided  by  the  continuity  equation, 
which  reduces  to  =  0  at  y— 0.  The  solution  for  v  and  u>y  together  with  the  con¬ 
tinuity  equation  and  the  definition  of  u>y  provide  equations  for  u  and  w.  Since  the 
pressure  is  not  computed  as  part  of  the  solution  algorithm,  it  is  computed  separately 
using  a  Poisson-solver  when  needed. 

The  turbulence  statistics  for  the  2D  data  are  in  agreement  with  existing  experi¬ 
mental  results,  and  have  been  widely  used  in  turbulence  research  (Moin,  1987).  It  is 
therefore  a  well  verified  code. 

2.3  Parameters  and  grid 

The  effect  of  periodic  boundary  conditions  on  turbulence  statistics  is  believed  to 
be  minimal  if  the  horizontal  dimensions  of  the  computational  box  are  sufficiently  large 
that  all  correlations  decay  t,o  zero  in  a  distance  less  than  half  the  computational  box 
size  in  each  of  the  horizontal  directions.  In  the  present  study,  the  physical  dimensions 
of  the  computational  domain  are  ( Lx ,  Ly,  Lz)  ~  (4tt6,  26,  j7t6),  where  8  is  the  channel 
half  width.  The  number  of  grid  points  is  (NX,  NY,  NZ)  =  (128,129,128).  The 
Reynolds  number,  based  on  6  and  the  2D  wall  shear  velocity,  ur2D  =  ‘s  180, 
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which  corresponds  to  a  Reynolds  number  based  on  the  centerline  velocity  of  3300.  For 
the  2D  flow,  the  grid  dimensions,  in  wall  units,  are  Ax+  =  17.67  and  A z+  =  5.89, 
and  the  first  grid  point  off  the  wall  is  at  y+  —  0.054. 

The  governing  equations  are  normalized  with  the  2D  shear  velocity,  uT 2D,  and  6, 
which  gives  an  initial  non-dimensional  streamwise  pressure  gradient  ^  =  —  1.  For 
the  3D  simulation,  ^  is  set  to  keep  the  mean  mass  flux  in  the  x-direction  constant. 
The  imposed  spanwise  pressure  gradient  is  set  to  —10.  For  the  time-period  considered 
in  this  study,  the  mean  flow  angle  was  less  than  45°,  which  is  low  enough  for  the  3D 
flow  to  be  sufficiently  resolved  with  the  2D  grid.  Some  simulations  were  performed 

q  D 

with  an  extremely  high  spanwise  pressure  gradient,  =  —100.  This  was  done  in 
order  to  make  the  effects  large  and  thus  easier  to  study.  For  these  simulations  it 
was  necessary  to  use  NX= 256  because  of  the  large  turning  angle  which  required  the 
horizontal  grid  to  be  closer  to  a  square. 

2.4  The  minimal  flow  unit 

One  of  the  difficulties  with  the  study  of  organized  structures  in  turbulent  bound¬ 
ary  layers  is  the  presence  of  several  structures  located  randomly  in  space  and  time. 
The  minimal  flow  unit  (Jimenez  and  Moin,  1990)  is  a  channel  flow  in  which  the 
dynamics  is  forced  to  be  much  simpler.  The  horizontal  dimensions  in  the  channel 
flow  code  are  reduced  to  the  smallest  size  that  can  maintain  turbulence.  Only  one 
low-speed  streak  and  often  just  two  or  three  vortical  structures  are  found  within  the 
computational  domain.  The  flow  is  called  a  minimal  unit  because,  apparently  with 
a  minimum  set  of  structures  the  basic  cycle  of  turbulence  generation  and  dissipation 
is  maintained.  Even  though  the  flow  would  not  occur  naturally,  it  does  satisfy  the 
Navier-Stokes  equations  and  the  turbulence  statistics  are  very  similar  to  those  found 
in  a  “regular- sized”  channel.  It  is  therefore  believed  that  the  key  mechanisms  that 
maintain  turbulence  in  the  minimal  flow  are  nearly  the  same  as  that  in  the  regular 
channel. 

In  the  minimal  channel,  the  spanwise  dimension  of  the  computational  box  was 
slightly  less  than  100  wall  units  and  its  streamwise  dimension  was  about  250-350  wall 
units.  These  correspond  approximately  to  the  mean  spanwise  distance  betw’een  the 
streaks  and  the  streamwise  extent  of  vortical  structures  in  the  channel  flow. 
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The  minimal  flow  studied  in  this  work  had  a  Reynolds  number,  Re i  —  = 

2000  and  physical  dimensions  (Lz,  Ly,  Lz)=(t8,  28,  ^8)  while  the  number  of  grid- 
points  were  (NX,  NY,  NZ)— (32,129,16),  where  Ui  is  the  laminar  center  line  velocity 
with  the  same  mass  flux. 

In  the  study  to  be  described  in  §3.1.3,  the  minimal  channel  was  used  primarily 
to  obtain  an  increased  understanding  of  the  dynamics  of  the  “2 D”  wall  turbulence. 
In  addition  to  the  advantage  of  a  simplified  turbulence,  the  small  size  of  the  com¬ 
putational  domain  makes  it  much  easier  to  follow  a  single  structure  in  time.  Due 
to  periodic  boundary  conditions,  the  structure  of  interest  will  remain  throughout  its 
lifetime  inside  the  computational  domain  which  is  so  small  that  it  can  easily  be  found. 

A  spanwise  pressure  gradient  was  applied  to  the  minimal  flow  unit  as  well,  but 
because  the  spanwise  size  of  the  computational  domain  is  very  small,  the  comput  tions 
were  meaningful  only  for  very  small  turning  angles.  The  spanwise  pressure  gradient 
was  therefore  chosen  to  be  lower  than  in  the  regular  3D  channel,  l^rl  =  5|^p|.  It 
was  also  useful  to  apply  spanwise  pressure  gradients  of  both  signs  to  the  same  initial 
field  to  see  how  the  same  eddy  would  respond  to  mean  streamwise  vorticity  of  both 
signs. 


2.5  Preliminary  considerations  for  the  3D  turbulent  channel  flow 


The  purpose  of  this  section  is  to  show  that  the  3D  channel  flow  considered  in  this 
study  behaves  like  a  shear-driven  TBL  and  to  describe  its  global  characteristics.  It  is 
useful  to  first  look  at  the  laminar  case.  The  governing  equations  are  then 


duj  dp  1  d2Ui 

Hi  ~  ~dHi  +  TteUyt 


(2.5.1) 


The  equations  are  uncoupled  and  the  uj -profile  will  not  change  when  the  spanwise 
pressure  gradient  is  imposed.  Because  w  is  independent  of  u ,  it  will  evolve  as  an 
impulsively  started  2D  channel  flow.  It  is  more  instructive  to  consider  the  governing 
equation  for  wy  =  ^  and  its  boundary  conditions. 


dwy  1  d2Wy 

HT  =  Yedy2 


(2.5.2) 


dwy  _  ?P 

dy  y=0  dz 


=  0. 
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where  y  —  6  is  the  channel  centerline.  It  can  be  seen  that  the  spanwise  pressure 
gradient  decs  not  influence  ^  inside  the  flow;  it  generates  a  flux  at  the  wall,  which 
diffuses  outward.  A  spanwise  boundary  layer  will  therefore  grow  outward  from  the 
wall.  Outside  this  boundary  layer,  w  has  no  spatial  gradients,  and 


woo (t)  ~  ~  j 


dP 

dz 


0 


A  series  solution  for  w  can  be  found  for  the  entire  domain  using  the  method  of 
separation  of  variables 

dp  1  o  l(f+n)Tl  4  1  — 3  1 

w(t,y)  = -~  Re  {y  -  ~y2  -  ,[(-  +  n)jr]  sin[(- +  n)jry]}  (2.5.3) 


n=0 


A  similarity  solution  exists  for  the  case  of  infinite  distance  between  the  walls 

w{y)  ^ t  {c  (r)2  +  \)  f  {y2  +  \)  7t~^dT]  -rj2}  (2.5.4) 


where, 


CO 


1  =  J  (y2  +  1)  2e  drj ,  rj  =  \J Re/2~^- 

0 


This  solution  is  in  excellent  agreement  with  eqn.  (2.5.2)  as  long  as  the  spanwise 
boundary  layer  is  thin  compared  to  the  channel  half  width. 

When  the  flow  is  turbulent,  the  governing  equations  for  the  mean  velocities,  (7,. 
are 

dUj 
dt 


dP  1  d2Ui  du'-v’ 
-  +  *  * 


(2.5.5) 


dii  Re  dy ^  dy 
As  in  the  laminar  flow,  the  vertical  mean  velocity  component  is  zero.  It  can  be  shown 
that  the  mean  pressure  is 


P(x,y,  z,t)  ~  ox  +  flz  -  pv°{yj) 


(2-5.6) 


plus  an  arbitrary  constant,  where  a  and  /?  are  the  imposed  horizontal  mean  pressure 
gradients.  Their  values  in  the  present  flow  were  given  in  §2.3.  Since  v'w'  is  zero  ini¬ 
tially,  W  will  grow  like  the  laminar  flow  until  v'w1  becomes  significant.  The  spanwise 
boundary  layer  will  grow  outward  from  the  wall  in  the  turbulent  case  as  well,  while 
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further  out  W(t)  =  independent  of  y.  The  mean  velocity  in  the  x  direction, 

U ,  is  affected  by  the  spanwise  pressure  gradient  through  only.  Any  change 

in  U  is  therefore  due  to  a  change  in  the  structure  of  the  turbulence.  When  the  flow 
relaxes  back  to  a  new  steady  state,  the  Reynolds  stresses  will  be  higher  than  that  of 
the  initial  flow  since  the  Reynolds  number  is  higher. 

The  pressure  driven  3D  channel  flow  considered  here  is  equivalent  to  a  channel 
flow  with  walls  moving  with  a  uniform  spanwise  acceleration.  This  can  be  shown  as 
follows.  The  transformation 


.  dP, 

W  —  W - r — t 

az 


„  ldP  2 
Z~Z  2  dz  1 


(2.5.7 a) 


dw  „  dw  _  dp'  1  d^w 
dt  dz  Re  dy^ 


(2.5.76) 


tv(y  =  0,  t)  =  w(y  =  26,  t)  =  —  t  (2.5.7c) 

This  is  the  governing  equation  and  boundary  conditions  for  a  channel  flow  whose 
walls  are  accelerating  in  the  spanwise  direction  at  a  constant  rate. 


2.6  3D  channel  flow  vs  3D  boundary  layers 


The  3D  channel  flow  studied  here  should  be  contrasted  to  3D  boundary  layers 
studied  experimentally.  As  mentioned  in  Chapter  1,  most  pressure-driven  spatially 
evolving  3D  boundary  layers  encounter  inviscid  skewing  of  the  vorticity  vector  in  the 
outer  flow.  This  effect  is  absent  in  the  present  simulations  where  all  3D  effects  diffuse 
outward  from  the  wall.  The  present  study  examines  the  effect  of  three-dimensionality 
on  the  near  wall  turbulence.  Most  3D  boundary  layers  studied  experimentally  or 
encountered  in  applications  are  statistically  stationary  whereas  the  3D  channel  flow 
of  the  present  study  is  temporally  evolving.  Moreover,  due  to  spatial  homogeneity,  at 
any  instant  the  ensemble  averaged  streamlines  are  straight.  However,  the  flow  features 
are  similar  if  the  streamwise  coordinate  and  time  are  related  via  Taylor’s  hypothesis. 
Figure  2.1  shows  the  mean  flow  particle  paths  in  the  channel  flow,  showing  the  typical 
curvature  and  the  wall-normal  skewing  expected  in  spatially  evolving  boundary  layers. 
Unfortunately,  in  many  experimental  studies  of  3D  boundary  layers,  the  means  for 
generating  the  three-dimensionality  has  also  introduced  adverse  streamwise  pressure 
gradients.  Some  of  the  observed  effects  may  be  attributed  to  the  adverse  pressure 
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gradient  rather  than  to  three-dimensionality.  The  3D  channel  flow  studied  here  does 
not  have  such  additional  complex  effects. 

2.7  Methods  of  identifying  the  mechanics  underlying  the  3D  effects 

In  this  section  various  statistical  methods  and  a  technique  for  the  analysis  of  the 
unsteady  flow  dy nannies  are  described. 

The  probability  density  function  (PDF)  is  a  useful  quantity  for  linking  changes 
in  statistical  quantities  to  changes  in  the  turbulence  structures.  PDFs  show  the 
distribution  of  a  quantity,  /,  and  are  normalized  such  that 

P(f)df  =  1.  (2.7.1) 

where  P(f)  is  the  PDF  of  /.  When  the  PDF  is  multiplied  (weighted)  by  some  function 
g{f),  the  integral  over  all  values  of  /  is  equal  to  g.  A  plot  of  the  weighted  PDFs  for 
the  2D  and  the  3D  flows  will  therefore  show  how  much  the  change  in  occurrence 
of  each  /-value  contributes  to  the  change  in  g.  The  joint  PDF  between  two  flow 
variables,  f\  and  /2,  shows  how  the  different  combinations  of  /i-values  and  /2-values 
are  distributed.  When  it  is  weighted  by  some  function  gifx^h),  it  shows  how  each 
combination  contributes  to  g.  In  the  current  work,  this  is  used  to  study  how  each 
combination  of  u \  and  u'-  contributes  to  the  change  in  the  Reynolds  shear  stress  uj-u'-. 

The  quadrant  analysis  (Willmarth  k  Lu  1972;  Wallace  et  al.  1972)  of  tt'-u'  is  a 
useful  supplement  to  weighted  joint  PDFs.  Each  value  of  u'-u'j  is  classified  into  one  of 
four  quadrants  according  to  the  sign  of  and  u'y  The  Reynolds  shear  stress  can  then 
be  separated  into  a  sum  of  four  terms,  each  equal  to  the  integral  of  u'uyP(uj-,  u'  )  over 
that  quadrant.  The  change  in  a  term  shows  the  integral  effect  of  all  changes  within 
the  quadrant  and  how  it  varies  with  the  distance  from  the  wall. 

Simpson  and  Devenport  (1990)  pointed  out  the  usefulness  of  an  octant  analysis 
in  the  study  of  3D  TBLs.  The  Reynolds  stresses  are  separated  into  8  terms  according 
to  the  sign  of  each  of  the  velocity  fluctuations.  For  example,  the  first  octant-term  for 
u'v1  is  equal  to  the  sum  of  v!v'  over  all  points  where  u',  v1  and  w'  are  all  positive, 
divided  by  the  total  number  of  points.  We  will  use  the  octant  analysis  to  distinguish 
between  the  contributions  to  the  Reynolds  stresses  by  streamwise  vortices  of  positive 
and  negative  signs  of  rotation.  When  the  octant-contribution  for  the  2D  flow  is 
subtracted,  it  can  be  seen  how  each  octant  contributes  to  the  change  in  the  Reynolds 
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stress.  In  the  current  work,  only  the  change  in  magnitude  of  the  contribution  for  each 
octant  is  considered. 

The  average  flow  field  around  an  event  (loosely  defined  as  a  fluid  motion  in 
some  specified  direction)  or  a  structure  can  be  computed  using  conditional  averaging 
(Blackwelder  &  Kaplan  1976;  Kim  k  Moin  1986).  We  have  examined  the  effect  of 
three-dimensionality  on  conditional  eddies.  In  the  current  work  this  is  done  in  the 
following  manner  :  A  range  is  specified  for  a  chosen  quantity  at  a  particular  distance 
from  the  wall.  For  each  point  which  satisfies  these  conditions,  the  origin  of  the  coor¬ 
dinate  system  is  moved  to  that  point  and  the  surrounding  flow  field  is  added  to  that 
from  the  previous  points. 

Two-point  correlations  are  used  to  extract  information  about  the  size  and  rela¬ 
tive  locations  of  structures.  The  direction  of  strongest  correlation  was  used  to  esti¬ 
mate  the  average  angle  of  structures  as  follows.  In  a  horizontal  plane,  the  angle  of 
maximum  correlation,  8r,  is  derived  from  £;f{x,z)f(x  -j-  Ai,:  +  ctAx)  =  0  where 
a  =  tan(0ft)  =  .  A  Taylor  expansion  and  using  the  horizontal  flow  homogeneity 

and  periodicity  of  this  flow  leads  to 


a/ a/  a/ 5/ 

dx  dz  dz  dz' 


(2.7.2) 


Finally,  a  non-statistical  approach  called  the  companion  flows  method  was  used 
to  study  the  unsteady  flow  dynamics.  Two  simulations,  one  with  and  one  without  a 
spanwise  pressure  gradient,  were  started  from  the  same  “2D”  initial  field.  Gradually, 
the  two  flows  deviate  due  to  3D  effects.  Comparison  of  the  two  flows  allow  identifica¬ 
tion  of  3D  effects  on  the  instantaneous  structures.  For  example,  the  flow  field  around 
a  vortex  in  the  “2D”  flow  can  be  compared  with  the  field  around  the  same  vortex  in 
the  3D  flow  and  those  vortices  which  generate  significantly  different  fluctuations  can 
be  detected.  Comparison  is  most  useful  for  early  times,  when  large  differences  are 
restricted  to  a  few  locations  in  the  flow. 


22 


CHAPTER  3 


NEAR- WALL  TURBULENCE  STRUCTURES  IN  “2D”  BOUNDARY  LAYERS 

We  refer  to  “organized  structures”  as  the  non-random  component  of  the  velocity 
fluctuations  in  a  turbulent  flow.  To  quote  Lumley,  an  organized  structure  is  “a 
pattern  which  is  recognizable  and  recurs  throughout  the  flow”.  They  are  known 
to  be  the  primary  contributors  to  the  Reynolds  shear  stresses  and  other  statistical 
quantities  and  are  also  believed  to  be  important  in  the  dynamics  which  maintains 
turbulence.  Several  types  of  structures  have  been  identified  in  boundary  layers,  but 
in  this  work  our  discussion  tends  to  focus  on  the  three-dimensional  effects  on  the 
widely  observed  quasi-streamwise  vortices,  high /low-speed  streaks  and  shear  layers. 
Streamwise  vortices  are  associated  with  large  values  of  v  and  w  and  u>*,  streaks  are 
associated  with  large  streamwise  velocity  fluctuations,  u',  normal  vorticity  u>y,  and 
spanwise  vorticity  fluctuations  cj'2;  whereas  shear  layers  protruding  from  the  wall  are 
associated  with  large  uz  away  from  the  wall.  A  brief  review  of  the  existing  knowledge 
and  new  findings  regarding  these  structures  are  given  below.  This  study  is  limited 
to  the  near-wall  region  because  the  spanwise  boundary  layer  originates  and  grows 
outward  from  the  wall. 

3.1  Quasi-streamwise  vortices 

A  vortex  is  a  swirling  motion  in  coordinates  moving  with  and  aligned  with  its 
axis.  Several  vortex  lines  coming  together  in  a  compact  region  of  space  is  also  a  good 
indicator  for  the  presence  of  a  vortex  in  that  region.  Both  definitions  will  be  used 
below.  The  velocity  vectors  in  any  vertical  cut  normal  to  the  mean  flow  direction  show 
that  quasi-streamwise  vortices  are  abundantly  present  in  the  near  wall  region  and  are 
associated  with  intense  velocity  fluctuations  (Figure  3.1a).  Streamwise  vortices  have 
been  shown  to  be  responsible  for  turbulence  generation  and  the  so  called  bursting 
event  is  merely  due  to  the  passage  of  such  vortices  past  the  measuring  station  (Kim 
&  Moin,  1986).  In  the  near-wall  region  (y+  <  20)  the  most  intense  v',  w’  and  u1  >  0 
are  associated  with  the  vortices  closest  to  the  wall  (Figure  3.1). 
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3.1.1  Average  size  and  location 

The  mean  location  of  the  near  wall  streamwise  vortices  is  at  y+  =  20,  (Kim, 
Moin  &  Moser,  1987).  Once  a  streamwise  vortex  was  identified  in  a  (z  —  y)  plane, 
its  streamwise  extent  was  determined  by  a  sequence  of  z  —  y  cuts  in  the  positive  and 
negative  x-direction.  It  was  determined  that  the  streamwise  extents  of  the  vortices 
rarely  exceed  Ax+  =  350.  This  is  consistent  with  a  mean  streamwise  spacing  between 
near- wall  vortical  structures  of  440  reported  by  Clark  and  Markland  (1971). 

3.1.2  Detection  of  vortices 


Both  definitions  of  a  vortex  in  terms  of  swirling  motion  and  collapsed  set  of 
vortex-lines  can  be  used  to  identify  vortices  in  instantaneous  flow  fields.  However, 
when  several  vortices  are  to  be  visualized  in  3D  space  it  is  useful  to  have  a  function 
whose  contour  surfaces  capture  vortices.  One  obvious  choice  for  such  a  function 
is  the  streamwise  vorticity.  A  study  of  how  well  contours  of  streamwise  vorticity 
fluctuations  capture  swirling  velocity  vectors  in  the  z  —  y  plane  reveals  that  contour 
levels  of  magnitude  60uT2d/6  does  a  fairly  good  job  (Figure  3.1a).  In  the  minimal 
channel,  streamwise  vorticity  fluctuations  of  magnitude  0.5 uiJS  were  used.  Low 
pressure  regions  was  used  by  Robinson  (1991)  to  detect  vortices,  but  in  the  near- wall 
region  of  the  channel  flow  no  threshold  on  pressure  was  found  which  captured  most  of 
the  vortices  present,  and  only  vortices.  The  pressure  minima  in  the  centers  of  near¬ 
wall  vortices  was  not  sufficiently  strong  compared  to  low-pressure  events  occurring 
very  close  to  the  wall.  One  way  to  compensate  for  this  is  to  identify  regions  of  local 
minima  for  p.  This  can  be  done  by  plotting  high  values  of  V2p.  This  procedure 
captured  vortices  slightly  better  than  contour  surfaces  of  ux.  Consider  the  Poisson 
equation  for  pressure, 


du^  duj 
dxj  dxi 


The  dominant  source  term  inside  streamwise  vortices  was  found  to  be  ,  which 
is  the  product  of  the  component  terms  in  u>x.  Hunt,  Wray  &  Moin  (1988)  suggested 
using  the  second  invariant  of  the  velocity  gradient  tensor  together  with  low  pressure 
levels  as  the  detection  criterion  for  vortices.  This  did  not  seem  to  improve  significantly 
on  the  detection  of  vortices  close  to  the  wall,  but  may  be  useful  further  away  from 
the  wall. 
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With  a  few  exceptions,  surfaces  of  streamwise  vorticity  fluctuations  will  be  used 
to  identify  vortices  in  the  present  study  because  tJx  is  easily  computed  and  gives  the 
sign  of  rotation  for  the  vortex  as  well.  For  quantitative  studies,  the  swirling  velocity 
vectors  in  a  vertical  plane  normal  to  the  vortex-axis  were  used  to  identify  vortices. 


3.1.3  Dynamics  of  quasi-streamwise  vortices  in  TBLs 


The  dynamics  of  quasi-streamwise  near-wall  vortices  was  studied  using  the 
minimal  channel  (see  §2.4).  It  is  convenient  to  study  the  dynamics  in  this  flow  because 
the  computational  domain  is  small  and  it  is  significantly  more  economical  to  follow 
the  temporal  evolution  of  the  flow  structures. 

Any  viable  description  of  the  flow  dynamics  is  considered  complete  if  it  is  cyclic. 
In  the  following,  the  process  leading  to  the  formation  of  streamwise  vortices  will  be 
described,  i.e.  the  sequence  of  events  leading  to  the  roll-up  of  a  layer  of  vorticity  into 
a  streamwise  oriented  vortex  rod  will  be  outlined. 

Figures  3.2  and  3.3  show  the  roll-up  of  a  layer  of  streamwise  vorticity  into  a 
streamwise  vortex.  Initially,  a  compact  streamwise  vortex  (shown  by  dashed  lines) 
and  two  layers  of  streamwise  vorticity  of  opposite  sign  are  discernible.  Of  interest 
is  the  vortex  layer  near  the  wall.  The  layer  is  seen  to  roll-up  and  form  a  compact 
streamwise  vortex.  A  perspective  view  of  the  first  roll-up  is  shown  in  Figure  3.4. 

The  sequence  of  cyclic  underlying  events  for  the  production  of  this  streamwise 
vortex  are  as  follows.  We  begin  with  an  inclined  streamwise  vortex  near  the  wall 
(corresponding  to  the  vortex  marked  by  the  dashed  lines  in  Figure  3.2).  Consider  the 
governing  equation  for  the  streamwise  vorticity 


Dojx 

~Df 


du  du  du  _9 
—  U)x  +  U>y  —  -f  V  V  Ux , 


(3.1) 


The  tilting  terms  on  the  right  hand  side  of  eqn.  (3.1)  are  simplified  as  follows 


du  du  dv  du  dw  du 

y  dy  +  z  dz  dx  dz  dx  dy 


(3.2) 


Instantaneous  contour  plots  of  all  the  terms  in  the  u >x  equation  clearly  revealed 
that  the  regions  of  large  values  of  the  second  term  in  the  right  hand  side  of  (3.2) 
coincide  with  the  layer  of  streamwise  vorticity  shown  in  Figure  3.2  (see  Figure  3.5). 
Therefore,  }s  the  dominant  contributor  to  the  formation  of  the  streamwise  layer 

of  vorticity  shown.  The  tilting  of  the  parent  streamwise  vortex  with  respect  to  the 
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wall  leads  to  significant  ^  (Figure  3.6a).  The  layer  is  stronger  on  the  downwash 
side  of  the  parent  vortex  due  to  higher  values  of  ^  encountered  there  (Figure  3.2). 

This  is  how  the  streamwise  vortex  layer  is  generated.  The  layer  is  unstable 
(Jimenez  &  Orlandi,  1992)  and  rolls  into  a  streamwise  vortex.  Examination  of  the 
terms  in  eqn.  (3.1)  reveals  that  once  the  vortex  is  formed,  the  stretching  term, 
becomes  dominant  and  intensifies  the  vortex  (Figure  3.7).  The  inclination  of  the 

ft 

newly  formed  vortex  generates  horizontal  velocity  derivative,  (Figure  3.66). 

Note  the  velocity  derivative  along  the  axis  of  the  inclined  vortex,  w/,  is 

duj  _  uffUj  dui 
dx  i  ~  ukuk  dxj 

Computation  of  each  term  showed  that  the  dominant  term  is 

u>%  du 
u:ku)k  dx 

It  is  interesting  that  even  though  the  vortex  is  situated  near  the  wall  the  term 
involving  ^  has  less  contribution  to  stretching  than  the  term  involving 

Note  that  according  to  the  process  described  above  only  tilted  streamwise  vortices 
are  capable  of  self-generation  and  self-induced  stretching.  It  is  well  known  that 
a  streamwise  vortex  near  a  wall  is  associated  with  a  concentration  of  vorticity  of 
opposite  sign  at  the  wall.  It  is  important  to  distinguish  between  this  vorticity  and 
the  vortex  layer  of  opposite  sign  described  here.  In  fact,  Orlandi  and  Jimenez  (private 
communication)  have  demonstrated  that  the  wall  streamwise  vorticity  does  not  play 
an  important  role  in  the  near  wall  dynamics. 

3.2  High-  and  low-speed  streaks 

Contours  of  u'  in  horizontal  planes  near  the  wall  are  highly  elongated  in  the 
streamwise  direction,  particularly  in  the  viscous  sublayer.  Low-speed  streaks  ( u '  <  0) 
are  longer  than  high-speed  streaks  (Robinson,  1991).  The  former  will  usually  have  a 
width  between  Az+  =  20  and  40  and  their  lengths  can  exceed  A:r+  =  2000.  Contours 
of  positive  u'  cover  wider  regions  and  become  rapidly  shorter  with  increasing  contour- 
levels.  Their  widths  vary  significantly,  but  are  mostly  between  Az+  =  40  and  80. 
The  average  distance  in  the  spanwise  direction  between  two  streaks  of  the  same  type 
(according  to  the  sign  of  u1)  is  well  established  as  being  close  to  Az+  =  100.  Streaks 
contain  most  of  the  intense  u1 ,  utz  and  near  the  wall. 


28 


The  significant  streamwise  extent  of  ^-contours  may  be  due  to  a  sequence  of 
vortices  following  each  other,  pumping  high-speed  fluid  toward  the  wall  and  low- 
speed  fluid  away  from  the  wall.  Plots  of  the  velocity  vectors  in  z  —  y  planes  along 
a  randomly  chosen  low-speed  streak  support  this  scenario  (a  few  planes  are  shown 
in  Figure  3.8).  Note  that  the  vortices  seen  near  the  streak  in  the  z  —  y  planes  in 
Figure  3.8  have  alternate  signs,  and  that  the  distance  between  the  vertical  planes  is 
slightly  more  than  350  wall  units  (see  §3.1.1).  This  picture  is  consistent  with  the 
dynamics  of  vortex-generation  as  described  in  §3.1.3.  There,  it  was  argued  that  near¬ 
wall  vortices  are  generated  below  tilted  near-wall  vortices  of  opposite  sign.  Since  the 
newly  generated  vortex  is  closer  to  the  wall,  it  will  have  a  lower  streamwise  convection 
velocity  and  will  therefore  end  up  behind  the  “parent-vortex”.  A  sequence  of  vortices 
with  alternate  signs  and  varying  distances  from  the  wall  can  thus  be  created. 

As  pointed  out  above,  contours  of  positive  u'  tend  to  be  shorter  than  those  of 
negative  u' .  Marking  of  an  instantaneous  flow  field  according  to  the  sign  of  v  reveals 
a  streaky  pattern  (Figure  3.9a).  However,  contours  become  rapidly  less  elongated 
when  the  magnitude  of  the  contour  level  is  increased  (not  shown)  and  both  high  speed 
and  low  speed  streaks  contain  regions  of  intense  v  which  are  not  highly  elongated. 
The  relationship  between  the  magnitude  of  u'  and  v'  can  be  seen  from  the  joint 
probability  density  function  between  u'  and  v'  (not  shown).  It  shows  that  in  the 
viscous  sublayer  the  magnitude  of  v!  >  0  is  correlated  with  the  magnitude  of  v'  <  0 
while  the  magnitude  of  u'  <  0  appears  to  be  independent  of  the  magnitude  of  v'  >  0. 
This  is  seen  even  more  clearly  when  the  joint  PDF  is  weighed  by  u/2  (Figure  3.10). 
Therefore,  high  speed  streaks  tend  to  be  shorter  than  low  speed  streaks  (Figure  3.96). 

The  difference  in  correlation  with  v'  between  u1  <  0  and  v!  >  0  is  explained  as 
follows:  A  high  magnitude  of  v'  <  0  near  the  wall  indicates  that  the  fluid  is  from 
a  region  far  away  from  the  wall  and  will  therefore  generate  high  u'  >  0.  On  the 
other  hand  regions  with  v'  >  0  pump  low-speed  fluid  from  a  rather  uniform  bed  of 
low-speed  fluid  (Figure  3.11). 

High-speed  streaks  tend  to  be  wider  than  low-speed  streaks  because  the  fluid 
moving  toward  the  wall  is  redirected  in  the  spanwise  direction  by  the  streamwise 
vortices  and  the  wall.  A  spanwise  movement  of  high-speed  fluid  will  generate  a  wider 
region  of  u'  >  0  while  a  spanwise  motion  of  low-speed  fluid  near  the  wall  does  not 
generate  u'  <  0  (outward  motion  does). 
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3.3  Near-wall  shear  layers 


When  the  vortices  eject  low-speed  fluid  away  from  the  wall,  a  strong  shear  layer 
develops  between  the  low-speed  fluid  and  the  high-speed  fluid  it  replaces.  Thin  regions 
of  ^  «  0  are  therefore  seen  on  top  of  low-speed  fluid  (Figure  3.12).  The  abundance 
of  shear  layers  protruding  from  the  wall  in  transverse  views  of  turbulent  boundary 
layers  was  reported  by  Jimenez,  Moin,  Moser  and  Keefe  (1988).  Johansson,  Her  and 
Haritonidis  (1987)  pointed  out  that  shear  layers  have  long  life  times.  Together  with 
the  high-speed  streaks,  the  shear-layers  contain  the  strongest  vorticity  fluctuations 
(u>'z)  in  a  TBL.  Robinson  (1991)  observed  that  these  shear-layers  could  “shed” 
spanwise  vortices.  This  is  not  believed  to  be  a  significant  event  for  the  near-wall 
turbulence  which  is  the  region  studied  in  this  work. 


(a) 


f 


Figure  3.1.  (a)  Sample  of  velocity  field  and  regions  of  \ux\  >  60 ur/6  (grey),  (b) 

regions  below  =  10  in  (a)  with  |u|  >  0.5ur,  (c)  |w|  >  ur  and  (d)  u'  >  2 ur  (black) 
and  u'  <  — 2ur  (grey).  The  Figure  demonstrates  that  vortices  closest  to  the  wall 
generate  the  most  intense  |u|,  |u;|  and  u'  >  0,  but  not  always  the  most  intense  u'  <  0. 
Tic  mark  are  spaced  0.25  apart. 
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Figure  3.2.  Time  evolution  of  contours  of  u)x  in  coordinates  moving  with  the 

structures.  Increment  is  0.25ti^/6  starting  with  0.5ul/6  (  -  )  and  -0.25 u^/S 

starting  with  -0.5 u^/S  ( - ).  Velocity  vectors  are  included  for  the  first  frame.  Tic 

marks  are  0.2  apart.  The  time  sequence  is  a,b,c,d,e,f 


Figure  3.4.  Time  evolution  of  contours  of  ux  —  0 .Su^/S  in  coordinates  moving  with 
the  structures.  The  z  —  y  cuts  across  the  contour  surfaces  are  also  seen  in  Figure  3.2. 
Tickmarks  are  0.2  apart.  The  time  sequence  is  a.,b,c,d,t,f. 
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Figure  3.5.  Components  o:  tilting  terms  for  the  streamwise  vorticity  associated  with 
the  vortex  formation  in  Figure  3.2.  Contour  levels  of  —  ( - )  and  cF 

( - )  start  at  0.05 (u/7^)2  and  increment  is  0.05(ujr,/^)2.  Contour  level  of  |u^x  I  is 

0.5 ui/S  (-  ■■  Tic  marks  are  0.2  apart. 
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Figure  3.6.  Horizontal  gradients  associated  with  a  tilted  near-wail  vortex.  Contours 

of  ux  =  Q.5ui/8  ( . - . . )  show  the  newly  formed  vortex  in  the  last  frame  in  Figure 

3.4.  Iso-contours  of  positive  ( - )  and  negative  ( - )  values  of  (a)  w  and  (6)  u 

are  inclined  due  to  the  tilting  of  the  vortex.  Tic  marks  are  spaced  0.25  apart. 


(a)  (b)  (c)  (d) 


Figure  3.8.  Velocity  vectors  in  vertical  planes  along  a  randomly  chosen  low-speed 
streak  at  y^D  =  5  (top,  black:*/  <  —  ur).  Vertical  lines  show  the  location  of  the  long 
low-speed  streak  in  the  top  figure.  The  horizontal  lines  is  the  plane  of  the  streak 
(y+  =  5).  This  illustrates  that  streaks  are  maintained  by  a  system  of  staggered 
vortices.  Tic  marks  are  0.25  apart  except  on  the  horizontal  axis  of  the  top  figure 
where  they  are  0.5  apart. 
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Figure  3.9.  (a)  Regions  of  positive  v  (grey)  at  y%D  =  5.4.  In  ( b )  grey  represent 

u'  >  0  and  the  lines  are  contours  of  u 1  =  ±I.5uT.  The  contours  of  u'  <  —  1.5«r  are 
longer  and  thinner  than  those  of  u!  >  1.5ur-  Tic  marks  are  0.5  apart. 
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Figure  3.11.  Schematics  showing  the  regions  of  origin  for  fluid  crossing  =  5 
over  a  time  A t.  The  region  of  fluid  swept  toward  the  wall  originate  from  regions 
extending  across  the  mean  velocity  gradient  while  the  ejected  fluid  originates  from  a 
region  extending  in  the  horizontal  directions. 
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layers  on  top  of  a  low  speed  streak.  (6)  lines  represent  u'  =  —  uT  and  the  grey  area  show  the  region  over  which  the  shear 
layers  are  accumulated  in  (a).  Tic  marks  are  0.5  apart  on  the  horizontal  axes  and  0.25  apart  on  the  vertical  axes. 


CHAPTER  4 


RESULTS  AND  DISCUSSION 

The  statistics  of  the  3D  flow  studied  here  were  reported  by  Moin  et  al.  (1990). 
They  found  a  reduction  in  the  total  Reynolds  shear  stress  (Figure  4.1),  turbulent 
kinetic  energy,  and  the  structure  parameter,  a  j.  The  reduction  in  turbulent  kinetic 
energy  was  due  to  a  large  decrease  in  u'2  and  a  smaller  decrease  in  v'2  (Figure  4.2a). 
In  coordinates  aligned  with  the  Reynolds  shear  stress  direction  (see  eqn.  2.16),  the 
spanwise  intensity  was  almost  unchanged  (Figure  4.26).  The  Reynolds  stress  budgets 
showed  that  a  reduction  in  the  production  terms  caused  the  drops  in  u/2  and  —  u'v' 
while  a  reduction  in  — caused  the  drop  in  v/2,  which  does  not  have  a  production 
term.  The  Reynolds  shear  stress  angle  lagged  behind  the  mean  velocity  gradient 
angle. 

In  this  chapter,  the  mechanics  underlying  these  changes  are  described.  Emphasis 
will  be  placed  on  the  evolution  of  the  six  Reynolds  stresses,  however,  evolutions  of 
the  mean  velocities  (§4.3),  vorticity  intensities  (§4.4)  and  the  two-point  correlations 
(§4.5)  are  also  discussed.  The  present  results  will  be  compared  with  the  experimental 
data  and  the  underlying  mechanisms  suggested  in  the  literature  (§4.6). 

It  is  known  (see  Chapter  3)  that  the  quasi-streamwise  vortices  are  the  primary 
contributors  to  the  Reynolds  shear  stress  in  2D  boundary  layers.  The  dynamics 
associated  with  these  vortices  will  therefore  be  studied  in  detail.  It  turns  out  that  the 
flow  dynamics  associated  with  positive  and  negative  streamwise  vortices  (the  same 
sign  as  or  opposite  to  the  mean  streamwise  vorticity)  are  different.  In  this  work 
the  sign  of  a  streamwise  vortex  will  be  defined  as  that  of  the  streamwise  vorticity 
fluctuation  in  the  center  of  the  vortex.  The  mean  x-vorticity  is  positive  in  the  present 
flow  due  to  the  direction  of  the  imposed  spanwise  pressure  gradient  (Figure  4.3). 

It  will  be  shown  below  that  the  key  mechanisms  for  the  reduction  in  the  Reynolds 
stresses  in  the  3D  flow  are  : 

I  Fluid  swept  toward  the  wall  by  positive  vortices  will  not  get  as  close  to  the  wall 
in  3D  as  in  2D,  and  will  therefore  generate  less  Reynolds  stresses.  In  coordinates 
moving  with  the  vortex  in  the  spanwise  direction,  the  fluid  will  be  forced  toward 
the  vortex-center  and  will  therefore  experience  a  lower  vortex-induced  vertical 
velocity  (Figure  4.4a). 


43 


II  Ejections  by  negative  vortices  generate  less  intense  velocity  fluctuations  because 
the  lifted  fluid  originates  at  a  higher  distance  from  the  wall  in  3D  (Figure 
4.46).  These  ejections  will  also  experience  higher  viscous  forces  leading  to  higher 
dissipation  in  the  turbulent  kinetic  energy  equation. 

These  are  the  only  significant  mechanisms  during  the  early  evolution  of  the  flow. 
Later  two  other  mechanisms  contribute  significantly  to  the  reduction  of  the  Reynolds 
stresses.  These  are: 

III  The  streamwise  vortices  are  shifted  with  respect  to  the  wall  layer  streaks  below 
them.  In  the  new  configuration,  the  vortices  pump  high-  rather  than  low-speed 
fluid  away  from  the  wall.  In  addition,  the  flow  trajectories  are  altered  with  adverse 
effect  on  the  Reynolds  shear  stress  (Figure  4.4c). 

IV  The  high  speed  fluid  swept  toward  the  wall  by  negative  vortices  is  convected  away 
from  the  vortex  in  the  spanwise  direction  by  the  spanwise  mean  flow  and  does 
not  reach  as  close  to  the  wall  as  in  2D  (Figure  4.4 d). 

Since  the  effect  of  the  spanwise  pressure  gradient  diffuses  outward  from  the  wall 
in  this  flow  (see  §2.5),  it  is  useful  to  study  the  initial  and  later  responses  separately. 
The  initial  period  is  defined  as  the  time  when  the  spanwise  boundary  layer  edge  is 
below  y%D  —  20  (y+  based  on  2D  wall  stress).  The  quasi-streamwise  vortices  are 
mostly  outside  the  spanwise  boundary  layer  and  the  flow  is  still  very  similar  to  a  2D 
channel  flow.  Comparison  of  the  2D  and  3D  flows  is  therefore  much  easier  than  at 
later  times  and  3D  effects  can  be  readily  identified.  One  of  the  difficulties  in  studying 
the  flow  during  the  later  period  is  that  the  there  is  no  obvious  choice  of  coordinates 
in  which  a  comparison  with  the  2D  flow  can  be  made  (see  §2.1). 

4.1  Reynolds  stresses  during  the  initial  period. 

For  some  time  after  is  imposed,  most  of  the  streamwise  vortices  will  be  outside 
the  spanwise  boundary  layer.  They  will  then  convect  in  the  spanwise  direction  without 
turning.  A  study  of  the  instantaneous  flows  showed  that  this  lasts  for  a  time  0.36 /uT2o 
when  =  — 10  (unless  otherwise  stated  all  quantities  are  normalized  with  the  2D 
wall  shear  velocity,  ut2d,  and  the  channel  half  width,  6).  This  is  the  time  it  takes  the 
spanwise  boundary  layer  to  reach  the  average  location  of  the  centers  of  the  near- wall 
vortices,  y^D  =  20.  The  lack  of  turning  at  f=0.3  can  be  seen  in  Figures  4.5  a  and  6, 
where  contours  of  u>'x  of  magnitude  00uT2D/6  (  which  in  §3.1.2  were  found  to  capture 
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vortices  fairly  well)  and  |^-|  >  50  at  the  wall  are  shown  for  the  two  companion 
flows  described  in  §2.7.  The  two  simulations  were  started  from  the  same  2D  initial 
field  t  =  0.3  earlier.  Both  vortices  and  streaks  have  been  convected  in  the  spanwise 
direction,  but  their  shapes  have  not  been  altered  significantly  by  the  spanwise  pressure 
gradient.  However,  the  turbulence  statistics  have  been  significantly  changed  within 
the  viscous  sublayer.  The  total  Reynolds  shear  stress  does  not  play  an  important 
role  in  the  viscous  sublayer,  but  a  drop  of  about  10%  in  this  region  indicates  that  a 
significant  change  in  turbulence  mechanics  has  already  occurred.  In  this  section,  the 
mechanics  causing  the  altered  statistics  at  2=0.3  will  be  discussed,  and  it  is  assumed 
that  the  mechanics  is  the  same  at  earlier  times. 

The  investigation  of  the  six  Reynolds  stresses  will  be  done  in  two  steps.  First, 
for  each  Reynolds  stress  the  events  which  are  contributing  most  to  its  change  are 
identified.  This  is  done  using  probability  density  functions  and  quadrant /octant- 
analysis  as  well  as  the  companion  flows  method  (§2.7).  Second,  these  events  are 
studied  in  detail  and  their  underlying  mechanics  are  identified. 

The  percentage  changes  in  the  Reynolds  stresses  are  used  to  identify  the  locations 
where  each  stress  is  most  affected.  At  2=0.3,  all  three  intensities  have  decreased  across 
the  3D  layer.  The  streamwise  intensity,  u'2,  has  dropped  by  2%  at  its  maximum 
location,  yfD  Rj  15,  and  the  maximum  percentage  drop  is  8%  at  y%D  «  4  (Figure  4.6). 
The  vertical  velocity  intensity,  o'2,  has  a  maximum  percentage  drop  of  9%  at  y?D  w  12 
while  the  spanwise  intensity,  u>'2,  has  a  maximum  reduction  of  9%  very  near  the  wall 
(Figures  4.7,  4.8).  The  streamwise  Reynolds  shear  stress,  — uV,  has  been  reduced 
by  5%  at  its  maximum  location  (y+p  «  30)  and  the  maximum  percentage  drop  is 
13%  very  close  to  the  wall  (Figure  4.9a).  The  spanwise  Reynolds  shear  stress,  v'w' . 
has  become  slightly  negative  with  its  maximum  magnitude  at  y%D  ss  10  and  the  ratio 
of  spanwise  to  streamwise  eddy- viscosity  is  0.3  across  the  viscous  sub-layer.  Finally, 
u’w'  has  become  slightly  positive  with  a  maximum  value  at  y%0  &  10  which  is  ten 
times  the  local  t/tu'-magnitude  and  half  the  local  uV-magnitude  (Figure  4.1a). 

4.1.1  Probability  density  functions,  octant-  and  quadrant  analysis 

The  intensity- weighted  PDFs,  tiJ2P(uj-),  show  which  ranges  of  velocity  fluctuation 
levels  are  contributing  most  to  the  aforementioned  reductions.  They  are  shown  in 
Figure  4.10  near  the  y-locations  of  maximum  percentage  changes,  but  the  trend  is 
virtually  the  same  across  the  viscous  sublayer.  The  streamwise  intensity  is  reduced 
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due  to  a  reduction  of  the  most  negative  and  medium-range  positive  u'  while  the 
vertical  intensity  is  reduced  due  to  reductions  of  medium-range  positive  and  negative 
v' .  Medium-range  positive  w'  and  a  wide  range  of  negative  values  causes  the  drop  in 

f\  / 

the  spanwise  intensity.  The  asymmetry  of  the  weighted  PDF  for  at  i= 0.3  may 
indicate  that  vortices  of  opposite  sign  are  affected  differently.  Note  that  the  PDF  of 
w'  very  near  the  wall  is  proportional  to  that  of  . 

According  to  the  quadrant  analysis  (Figure  4.11),  the  weakening  of  the  events  in 
the  second  quadrant,  u'  <  0  with  v'  >  0  and  in  the  fourth  quadrant,  u'  >  0  with 
v'  <  0,  contribute  about  equally  to  the  reduction  in  —u'v'. 

The  octant  analysis  (Figure  4.12)  gives  key  information  about  the  location  of 
weakened  events  (Simpson  &  Devenport  1990).  It  shows  that  in  the  viscous  sublayer 
the  reduced  quadrant-2  events  tend  to  be  correlated  with  positive  w'  while  the 
reduced  quadrant-4  events  tend  to  be  correlated  with  negative  w'.  Since  w'  is  mostly 
associated  with  streamwise  vortices  in  this  early  time-period,  this  indicates  that  those 
quadrant-2  events  that  are  weakened  are  generated  by  negative  vortices  while  reduced 
quadrant-4  events  tend  to  be  generated  by  positive  vortices  (see  Figure  4.3).  This  is 
in  accordance  with  mechanisms  I  and  II  presented  earlier  in  this  chapter.  Outside 
the  spanwise  TBL,  the  octants  with  u'v'  <  0  and  w'  >  0  are  most  reduced.  The 
octant  analysis  does  not  provide  enough  information  to  establish  the  sign  of  the 
vortices  which  generate  these  weakened  events,  since  in  this  region  these  events  may 
be  below  or  above  the  vortex-centers.  Figure  4.13  shows  that  the  same  two  octants 
are  responsible  for  the  reduction  in  the  velocity  intensities  as  well. 

The  weighted  joint  PDFs  between  u'  and  w'  at  y+D  —  10  (Figure  4.14)  show 
that  the  increase  in  u'w'  is  due  to  a  reduction  of  the  most  negative  u'  correlated 
with  positive  w' ,  and  of  medium-range  positive  u'  correlated  with  negative  w' .  This 
is  consistent  with  the  octant  analysis  for  u'v' .  The  quadrant  analysis  (Figure  4.15) 
shows  that  both  quadrants  which  generate  negative  u'w'  contribute  significantly  to 
the  increase  in  u'w' .  The  regions  of  reduced  positive  and  negative  u  are  therefore 
causing  the  change  in  this  Reynolds  stress  as  well. 

The  weighted  joint  PDFs  between  v '  and  w'  at  y$D  =  10  (Figure  4.16)  show  that 
the  negative  v'w'  is  due  to  a  reduction  of  intense  positive  and  negative  v'  correlated 
with  intense  w'  of  the  same  sign  (quadrant  1  and  3).  This  is  supported  by  the  quadrant 
analysis  (Figure  4.17).  The  reduced  values  of  w'  are  well  correlated  with  the  reduced 
values  of  v' .  Therefore,  negative  v'w'  is  generated  due  to  the  weakening  of  certain 
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events  which  is  very  different  from  the  mixing  mechanism  typically  advanced  for 
generation  of  negative  u'v’.  It  is  for  this  reason  that  eddy  viscosity  models  can  not 
predict,  at  least,  the  early  stages  of  three-dimensionality.  It  should  be  noted  that  the 
mechanisms  associated  with  events  I  and  II  (Figure  4.4)  described  earlier  have  now 
also  been  identified  with  generation  of  u'w1  and  —v'w'. 

4.1.2  The  companion  flows  method 

The  comparison  of  the  two  companion  flows  (see  §2.7)  will  now  be  used  to  confirm 
and  expand  these  findings.  Regions  with  ju'|  >  2,  which  according  to  the  PDF  of  u' 
are  most  affected  by  the  spanwise  pressure  gradient  at  ytD  =  5,  are  shown  for  the  two 
flows  in  Figure  4.18.  The  3D  flow  field  has  been  moved  in  the  negative  2-direction  to 
account  for  the  spanwise  convection  so  the  two  fields  become  easier  to  compare.  The 
shifted  streamwise  fluctuations  will  be  denoted  by  f/,  where  u'(x,z)  =  u'(x,z  +  /) 
and  /  =  l{y,t)  is  the  spanwise  length  that  the  3D  flow  has  been  convected.  /  was 
found  by  comparing  v!  in  the  two  instantaneous  flow  fields  at  one  j-locatior  (Figure 
4.19).  It  can  be  seen  that  the  two  flows  are  very  similar.  In  order  to  identify  the 
subtle  events  causing  the  changes  in  the  statistics,  the  regions  which  contribute  most 
to  the  change  in  i/2  are  located.  These  are  defined  as  the  regions  where  u'3D  —  u'2D  is 
most  negative.  In  Figure  4.20,  contours  of  the  most  negative  u3D  -  u'2D  at  y2D  =  5  are 
shown.  The  plot  also  shows  the  sign  of  w'2f>  (w'3D  has  mostly  the  same  sign  as  its  2D 
counterpart).  It  can  be  seen  that  the  regions  which  have  a  reduction  in  u'2  are  well 
correlated  with  w'  with  sign  opposite  to  u'  (mechanism  I  &  II).  A  corresponding  plot 
which  showed  the  sign  of  v'  instead  indicated  that  the  reduced  {/-regions  correlate 
well  with  v'  of  opposite  sign.  Figure  4.21  shows  the  regions  which  contribute  most 
to  the  reduction  of  each  of  the  velocity  intensities.  It  indicates  that  the  event’  which 
cause  the  reduction  in  u'2  are  also  causing  the  reduction  in  i>'2  and  te/2 . 

The  vortices  which  are  associated  with  reduced  ?/  can  be  seen  in  ;  —  y  plane;  cut 
through  the  regions  identified  in  Figure  4.21.  In  almost  every  case  the  reductions  in 
positive  u'  were  associated  with  posit-ve  vortices  and  the  reductions  in  negative  u' 
were  associated  with  negative  vortices.  The  strengths  and  locations  of  the  vortices 
varied,  but  most  were  located  between  y3D~ 20  and  40  and  were  of  low  or  medium- 
range  strengths.  A  sample  of  each  of  the  vortices  and  the  corresponding  patches  of 
the  affected  velocity  regions  are  shown  in  Figures  4  22  and  4.23  in  coordinates  moving 
with  the  vortices  in  the  spanwise  direction. 
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The  Reynolds  stresses  a.e  a'ered  at  <=0.3  because  near- wall  vortices  produce  less 
intense  velocity  fluctuations.  Wlir'  these  vortices  have  a  positive  si  ^  of  rotation, 
they  tend  to  generate  reduced  u'  >  0,  v'  <  0  and  w'  <  0  (sweep),  and  when  they 
have  a  negative  sign  they  tend  to  generate  reduced  u'  <  0,  v'  >  0  and  w'  >  0 
(ejection).  This  causes  the  velocity  intensities  and  ~u'v'  to  be  reduced  and  —  v'wf 
u'w'  to  become  positive  (§4.2.1).  The  mechanisms  underlying  these  changes  v  ’!!  now 
be  discussed. 

Particle  trajectories  about  negative  vortices 

When  fluid  particles  enter  the  viscous  sublayer,  their  paths  will  be  affected  by 
the  spanwise  mean  velocity  gradient.  This  change  in  trajectory  is  linked  to  the 
alteration  of  the  Reynolds  stresses.  In  what  follows  it  is  most  convenient  to  ->sider 
the  equivalent  flows  where  the  flow  becomes  3D  due  to  a  spanwise  acceleration  of  the 
walls  (eqn.  2.5.7),  because  then  the  vortices  will  not  move  significantly  in  the  spanwise 
direction.  “Streamlines”  within  the  z  —  y  plane  (jf  =  ^;)  WM  therefore  indicate  how 
the  particle  paths  are  affected  by  the  spanwise  TBL.  Figures  4.24  shows  “streamline” - 
segments  for  the  sample-vortices  in  Figures  4.22  and  4.23.  Lines  for  both  companion 
flows  are  shown  together,  traced  to  or  from  a  common  point.  The  discussion  below 
is  based  on  these  figures. 

In  2D,  vortices  with  a  negative  sign  of  rotation  will  move  fluid  particles  in  the 
positive  2-direction  below  their  centers.  In  3D  this  motion  is  opposed  by  the  spanwise 
mean  flow.  If  a  fluid  particle  is  sufficiently  close  to  the  vortex-center,  it  will  encounter 
a  small  spanwise  mean  motion  and  will  have  sufficient  momentum  to  move  in  the 
positive  2-direction  and  be  ejected  as  iu  2D.  However,  a  fluid  particle  which  in  2D 
moves  down  close  to  the  wall  and  is  then  ejected  (and  generates  u'  <<  0),  will  be 
redirected  by  the  moving  wall  and  dragged  away  from  the  vortex  along  the  wall  (see 
also  sketch  in  Figure  4.46).  In  3D  the  fluid  which  is  ejected  by  the  vortex  is  convected 
toward  the  vortex  by  the  moving  wall  and  will  tend  to  originate  further  away  from 
the  wall  than  the  fluid  ejected  in  2D  (Figure  4.246).  Recall  that  the  PDF’s  of  u1 
indicated  that  the  most  negative  values  of  u'  were  most  affected  in  3D  Such  large 
negative  values  are  generated  by  fluid  originating  from  the  closest  proximity  of  the 
wall.  In  Figure  4.246  the  fluid  particles  originating  below  the  vortex  and  closest  to 
the  wall  are  swept  to  the  left  in  3D  and  are  not  shown.  Only  “streamlines”  near 
the  region  of  significant  reduction  in  the  magnitude  of  ?/  are  shown.  The  strongest 
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vortices  close  to  the  wall  at  <=0.3  are  able  to  move  fluid  in  the  positive  2-direction 
below  y+D  =  5  as  well  as  lifting  fluid  convected  toward  them.  They  will  therefore 
generate  intense  u!  <  0  in  3D  as  in  2D.  In  2D  these  streamlines  originate  from  below 
the  vortex  and  in  3D  they  originate  from  the  far  right  of  the  vortex. 

Moin  et  al.  (1990)  reported  a  higher  dissipation  of  turbulent  kinetic  energy  in 
3D.  This  appears  to  be  due  to  the  spanwise  convection  of  the  ejected  low-speed  fluid 
over  high-speed  fluid  regions  and  thereby  increasing  the  magnitude  of  <  0  (Figure 

4.25) .  As  will  be  shown  later,  the  probability  density  function  of  vorticity  fluctuations 

o  / 

show  that  the  probability  of  large  negative  increases  significantly. 

Particle  trajectories  about  positive  vortices 

Vortices  with  a  positive  sign  of  rotation  will  in  2D  move  fluid  toward  the  wall 
in  the  negative  z-direction  below  their  centers.  In  3D  the  spanwise  TBL  will  drag 
this  fluid  toward  the  vortex-center,  where  the  vortex-induced  vertical  motion  is  lower. 
Therefore,  it  will  not  reach  as  close  to  the  wall  as  in  2D  (Figure  4.24a,  4.4a).  This 
reduces  the  magnitudes  of  both  positive  u'  and  negative  v' .  The  large  values  of  w' 
near  the  wall  are  due  to  redirection  of  the  fluid  swept  toward  the  wall.  Therefore,  the 
reduction  of  negative  w'  close  to  the  wall  is  consistent  with  this  scenario. 

This  concludes  the  dynamics  of  near-wall  fluid  motions  associated  with  positive 
and  negative  vortices.  In  addition,  as  seen  in  Figure  4.9  the  —  u'v'  is  also  reduced 
outside  the  spanwise  boundary  layer.  The  same  approach  was  followed  in  order  to 
identify  the  regions  of  reduced  —  u'v'  in  the  outermost  spanwise  boundary  layer.  At 
y+o  —  25,  the  regions  where  —u'v'  were  most  reduced  were  found  to  correlate  well  with 
w  >  0  which  is  consistent  with  the  octant  analysis.  Vertical  cuts  through  regions  of 
reduced  quadrant-2  events  show  that  they  are  associated  with  negative  vortices  with 
centers  further  away  from  the  wall  (such  that  w'  >  0  below  them,  see  sketch  in  Figure 

4.26) .  In  the  same  way  it  was  found  that  the  reduced  quadrant-4  events  at  y*D  =  25 
were  associated  with  positive  vortices  with  centers  closer  to  the  wall  (w1  >  0)(see 
sketch  in  Figure  4.26  &  Figure  4.27).  The  signs  of  the  vortices  which  are  associated 
with  reductions  in  these  quadrants  are  therefore  the  same  across  the  TBL.  That  is, 
negative  vortices  are  associated  with  quadrant  2  reductions  and  positive  vortices  with 
quadrant  4  reductions.  The  vortices  linked  to  the  reduced  quadrant-4  events  tend  to 
be  closer  to  the  wall  than  those  linked  to  reduced  quadrant-2  events  (see  sketch  in 
Figure  4.26). 
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4.2  Reynolds  stresses  during  the  later  period  of  three-dimensionality 

After  2  «  0.3,  the  flow  becomes  more  obviously  three-dimensional  in  the  sense 
that  the  streamwise  vortices  start  to  turn  in  the  direction  of  the  mean  spanwise 
pressure  gradient.  The  vortices  closest  to  the  wall  become  submerged  in  the  spanwise 
boundary  layer.  Figure  4.28  shows  surfaces  of  streamwise  vorticity  fluctuations  equal 
to  60uT2d/£  (which  in  §3.1.2  were  found  to  adequately  capture  streamwise  vortices) 
at  /=0.45  and  2=0.60.  Contours  of  |^-|  >  50  at  the  wall  are  also  included  to  show  the 
streaks.  The  turning  of  vortices  makes  the  flow  more  difficult  to  analyze,  because  the 
vortices  now  contribute  to  all  the  Reynolds  stresses.  In  2D,  u1  is  mostly  generated  due 
to  transport  of  high  or  low  speed  fluid  in  the  vertical  direction  by  the  vortices  whereas 
v'  and  w'  are  mostly  due  to  the  swirling  motion.  Because  the  vortices  have  varying 
turning  angles  according  to  their  size  and  distance  from  the  wall,  it  is  impossible  to 
find  a  unique  coordinate  system  for  comparison  of  the  2D  and  3D  flows  using  the 
companion  flows  technique. 

4.2.1  Changes  in  the  Reynolds  stresses 

At  2=0.9,  the  velocity  intensity  in  the  initial  mean  flow  direction,  u/2,  reaches 
its  maximum  percentage  drop  of  28%  at  y$D  «  8  and  recovery  starts  at  the  wall 
(Figure  4.6a).  When  the  coordinates  are  rotated  to  keep  the  streamwise  intensity  at 
a  maximum  (eqn.  2. Id),  u*2  starts  to  recover  at  2=0.6  and  the  maximum  percentage 
reduction  is  19%  (Figure  4.66).  The  vertical  intensity,  i/2,  starts  a  slow  recovery 
near  the  wall  soon  after  2=0.3  (Figure  4.7).  At  2=0.9,  the  vertical  intensity  has 
almost  recovered  to  its  initial  values  near  the  wall,  while  at  y+D  40  it  reaches  a 
maximum  percentage  reduction  of  19%.  The  spanwise  intensity,  u?'2,  increases  rapidly 
after  2=0.3  (Figure  4.8).  In  coordinates  rotated  to  keep  the  spanwise  intensity  at  a 
minimum  (eqn.  2. Id),  w remains  slightly  lower  than  the  2D  profile  in  the  time-period 
studied  (Figure  4.26).  The  turbulent  kinetic  energy  reaches  a  maximum  percentage 
drop  of  17%  at  2=0.7  and  ^  10  and  then  starts  to  recover.  The  Reynolds  shear 

stress,  — uV,  attains  its  maximum  percentage  drop  of  23%  at  y+D  %  30  and  2=0.9 

j - 2  - 2 

(Figure  4.6a).  The  total  Reynolds  shear  stress,  yuV  +  v'w'  , which  is  invariant 
with  respect  to  coordinate  rotation,  reaches  a  maximum  percentage  drop  of  21%  at 
2  =  0.8  (Figure  4.96).  Figure  4.1a  shows  the  increase  in  the  spanwise  Reynolds  shear 
stresses  —v'w’  and  u’w’  due  to  the  turning  of  the  structures,  u’w'  should  be  viewed  as 
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a  velocity  correlation  rather  than  a  shear  stress  because  it  does  not  contribute  to  the 
mean  momentum  equations.  The  Reynolds  shear  stress  angle  where  vrw't  —  0,  eqn. 
(2.16),  and  the  intensity  angle,  where  u*u4  =  0,  eqn.  (2. Id),  are  shown  in  Figure 
4.29.  It  can  be  seen  that  the  latter  angle  lags  behind  the  former. 

The  changes  in  the  turbulence  statistics  can  be  attributed  to  the  structural 
changes,  to  the  rotation  of  the  coordinate  system  aligned  with  the  streamwise  vortices, 
and  to  the  increased  mass  flow  rate.  Coordinate  rotation  effects  can  be  eliminated  by 
considering  quantities  which  are  invariant  with  respect  to  coordinate  rotation.  One 
may  also  consider  a  coordinate  system  which  minimizes  the  turning  effect  on  the 
flow  quantity  being  examined.  There  are  two  characteristic  angles  for  the  Reynolds 
stresses:  the  Reynolds  shear  stress  angle,  eqn.  (2.16),  and  the  intensity  angle,  eqn. 

(2,1  d).  When  the  coordinates  are  rotated  with  the  stress  angle,  we  have  v'w^=  0 
-  /  2  2 

and  =  y  u'v1  +v'w'  .  A  'Station  of  the  coordinates  with  the  intensity  angle 

maintains  u^F  at  a  maximum,  u/*2  at  a  minimum  and  u*u4=0.  In  the  latter  coordinate 
system,  for  example,  it  would  be  legitimate  to  compare  the  2D  and  3D  PDFs  of 
streamwise  velocity  fluctuations.  The  increased  mass  flow  leads  to  increased  Reynolds 
stresses  and  flow  recovery. 

In  the  following  section  the  structural  changes  causing  the  reduction  in  the 
Reynolds  stresses  will  be  identified  and  the  underlying  mechanisms  will  be  discussed. 
Most  of  the  studies  were  conducted  at  t=0.6  when  many  statistical  quantities  were 
near  their  maximum  change. 

4.2.2  Identification  of  the  most-affected  events  in  the  late  period 

It  will  be  shown  in  this  section  that  the  reduction  of  the  Reynolds  stresses  is 
due  to  a  weakening  of  sweeps  and  ejections  by  vortices  of  both  signs,  but  as  in  the 
initial  period,  sweeps  by  positive  vortices  and  ejections  by  negative  vortices  are  most 
affected. 

In  coordinates  aligned  with  the  maximum  streamwise  intensity  (eqn.  2. Id)  the 
intensity-weighted  PDFs  (Figure  4.10)  show  the  same  trend  as  in  the  initial  period. 
That  is,  the  streamwise  intensity  is  reduced  because  there  are  fewer  medium-range 
positive  and  intense  negative  ti*,  the  normal  intensity  is  reduced  because  there  are 
fewer  medium-range  v'  of  either  sign,  and  a  reduction  of  negative  causes  the 
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reduction  in  spanwise  intensity.  The  latter  is  consistent  with  a  reduction  of  sweeps 
associated  with  positive  vortices. 

The  quadrant  analysis  for  u*v'  at  *=0.6  indicates  that,  as  in  the  early  period 
(Figure  4.11),  both  sweeps  and  ejections  are  weakened  (Figure  4.30).  This  result  is 
independent  of  the  choice  of  coordinates.  In  general,  the  analysis  of  Reynolds  stresses 
containing  u't  and/or  v'  was  qualitatively  independent  of  the  choice  of  coordinates, 
while  stresses  containing  w+  were  strongly  coordinate-dependent. 

Conditional  averaging  (see  §2.7)  was  used  in  order  to  establish  how  the  changes  in 
sweeps  and  ejections  by  vortices  of  each  sign  contribute  to  the  changes  in  the  Reynolds 
stresses.  The  primary  difficulty  with  conditional  averaging  is  to  find  the  condition 
which  captures  the  structure  of  interest.  The  changes  at  *=0.3  can  be  used  to  test 
the  choice  of  condition,  since  the  most  affected  events  at  this  time  are  known  to  be 
sweeps  by  positive  vortices  and  ejections  by  negative  vortices.  Conditioning  on  QjL 
at  the  wall  has  the  advantage  that  no  subjective  choice  of  location  inside  the  flow  is 
needed.  The  range  20  <  <  50  and  —50  <  <  —20  gave  the  correct  trends 

(Figure  4.31).  A  narrower  range  (30-40)  also  revealed  the  correct  features.  In  the 
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later  period,  the  conditional  eddy  corresponding  to  =  ±(20  -  50)  at  the  wall  was 
used,  where  w J,  is  the  spanwise  velocity  in  coordinates  rotated  with  the  intensity  angle 
(eqn.  2. Id).  In  these  coordinates,  the  rms  of  is  at  a  minimum.  The  conditionally 
averaged  flow  shows  a  reduction  of  sweeps  and  ejections  by  vortices  of  both  signs,  but 
the  two  events  which  were  most  affected  during  the  initial  period,  sweeps  by  positive 
vortices  and  ejections  by  negative  vortices,  are  also  most  affected  during  the  later 
period. 

In  order  to  obtain  less  coordinate-sensitive  conditionally-averaged  data,  we 
experimented  with  conditioning  on  v.  If  a  vertical  motion  is  part  of  a  vortex,  it 
will  be  accompanied  by  a  vertical  motion  of  opposite  sign  either  to  its  right  or  to 
its  left  according  to  the  position  of  the  vortex.  Several  locations  and  ranges  for 
v  were  tested.  Unfortunately,  the  results  were  highly  dependent  on  the  choice  of 
conditioning,  but  those  which  gave  the  correct  trend  for  the  initial  period  showed 
the  same  trend  in  the  later  period  as  well.  Figure  4.32  shows  the  averaged  field  of 
negative  u*  around  —0.25  <  v  <  —0.2  at  y%D  =  10.  It  can  be  seen  that  the  negative 
values  of  u\  to  the  right  (positive  z)  are  most  reduced,  which  implies  that  negative 
vortices  are  most  affected  (Figure  4.3).  Note  that  the  u\  fields  corresponding  to  both 
positive  and  negative  vortices  are  depicted  in  the  conditionally  averaged  field. 
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The  octant  analysis  was  very  useful  in  the  initial  period.  In  the  later  period,  the 
coordinate-dependence  of  w *  poses  a  problem.  Although  a  study  of  the  instantaneous 
flow  fields  indicated  that  w+  in  coordinates  aligned  with  either  the  Reynolds  shear 
stress  (eqn.  2.1  b)  or  the  direction  of  maximum  intensity  (eqn.  2.1  d)  was  a  good 
indicator  of  the  sign  of  the  vortices,  the  respective  octant  analyses  led  to  qualitatively 
and  quantitatively  different  conclusions.  The  only  consistent  trend  resulting  from 
both  octant  analyses  was  strong  reductions  from  ejections  generated  by  negative 
vortices.  In  general,  there  were  reductions  of  sweeps  and  ejections  by  vortices  of 
both  signs,  but  the  magnitudes  of  the  reductions  varied  strongly.  The  octant  analysis 
for  u+v1  in  coordinates  aligned  with  the  shear  stress  direction  is  shown  in  Figure  4.33. 
In  order  to  develop  a  less  coordinate-dependent  octant  analysis,  was  replaced 
by  the  streamwise  vorticity  fluctuation,  as  a  better  indicator  of  the  sign  of  the 
vortices.  However,  this  approach  was  unsuccessful  because  it  did  not  correctly  predict 
the  well-established  trends  for  the  early  period. 

In  conclusion,  sweeps  and  ejections  by  vortices  of  either  sign  generate  reduced 
fluctuations  in  the  later  period.  The  two  events  which  were  most  affected  during  the 
initial  period  (I  and  II  in  Figure  4.4)  continue  to  be  most  affected  in  the  later  period. 
In  addition  the  other  sweep  and  ejection  events  by  vortices  of  both  signs  are  weakened 
(III  and  IV  in  Figure  4.4). 

4.2.3  Underlying  mechanisms  in  the  late  period 

q  ; 

Conditional  averaging  (Figure  4.31)  and  the  intensity- weighted  PDF’s  for  at 
the  wall  (Figure  4.10c)  indicate  that  sweep-motions  by  positive  vortices  are  reduced 
in  the  later  period  as  well.  The  mechanism  is  believed  to  be  the  same:  in  3D,  fluid  will 
be  moved  toward  the  vortex-center  (in  coordinates  moving  with  the  vortex)  where 
the  vertical  velocity  is  lower,  and  it  will  therefore  not  reach  as  close  to  the  wall  as  in 
2D  (Figure  4.4a). 

The  weakening  of  ejections  associated  with  negative  vortices  are  believed  to  be 
caused  by  the  same  mechanisms  as  in  the  initial  period:  that  is,  the  fluid  closest  to  the 
wall  is  convected  in  the  spanwise  direction  instead  of  being  lifted  (Figure  4.46),  and 
the  ejected  fluid  experiences  more  viscous  dissipation  due  to  convection  of  low-speed 
fluid  on  top  of  high-speed  fluid  in  3D  (Figure  4.25).  This  is  supported  by  the  budgets 
for  the  Reynolds  stresses  which  show  continued  reductions  due  to  less  production  and 
more  dissipation  (Moin  et  al.  1990).  It  is  also  supported  by  the  weighted  PDF  for  u't 
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(Figure  4.10a)  which  shows  reductions  in  the  magnitudes  of  the  most  negative  values 
of  u*  which  result  from  the  transport  of  fluid  from  near  the  wall. 

Ejections  due  to  positive  vortices  generates  less  intense  fluctuations  in  the  later 
period.  Two  mechanisms  are  believed  to  be  responsible  for  this.  First,  the  low-  and 
high-speed  streaks  are  shifted  with  respect  to  the  vortices  above  them  (Figure  4.4. c). 
The  vortices  then  pump  higher-speed  fluid  away  from  the  wall  and  thereby  generate 
less  Reynolds  shear  stress  than  if  they  had  pumped  low-speed  fluid.  It  can  be  seen 
in  Figure  4.34c  that  in  3D  high-speed  fluid  indeed  resides  under  the  ejection  side  of 
the  vortex  shown.  A  second  mechanism  is  the  change  in  particle  trajectories.  Fluid 
which  is  ejected  will  simultaneously  be  convected  away  from  the  vortex  into  a  region 
of  weak  vortex-induced  vertical  motion.  This  is  seen  in  Figure  4.34a  where  “stream¬ 
lines”  within  a  vertical  plane  normal  to  a  vortex  indicate  the  trajectory  of  ejected 
fluid. 

Sweep-motion  by  negative  vortices  generate  lower  levels  of  velocity  fluctuations  in 
the  later  period  (Figure  4.4 d).  This  is  contrary  to  the  initial  period,  when  this  type 
of  event  generated  slightly  more  intense  fluctuations  than  in  2D  (see  octant  analysis 
in  Figure  4.12).  When  the  vortices  are  submerged  in  the  spanwise  boundary  layer, 
the  fluid  swept  toward  the  wall  will  simultaneously  be  swept  away  from  the  vortex 
and  will  not  reach  as  close  to  the  wall.  This  is  seen  in  Figure  4.35  which  shows 
selected  “stream-lines”  around  a  negative  vortex  in  the  3D  flow.  This  completes  the 
presentation  of  the  evidence  and  the  corresponding  mechanisms  for  the  affected  events 
sketched  in  Figure  4.4. 

There  are  two  characteristic  directions  for  the  Reynolds  stresses.  These  are  the 
direction  of  maximum  st.reamwise  intensity  (eqn.  2.1  d)  and  the  direction  of  maximum 
streamwise  Reynolds  shear  stress  (eqn.  2.16).  The  latter  leads  the  former  (Figure 
4.29)  and  the  reason  for  this  will  now  be  discussed.  The  difference  between  these 
angles  is  believed  to  be  due  to  different  turning  angles  of  the  vortices  according  to 
their  distance  from  the  wall.  The  vortices  which  are  closest  to  the  wall  generate 
the  highest  magnitudes  of  positive  u'  and  negative  v'.  Their  high  turning  angles 
are  therefore  responsible  for  high  stress-  and  intensity  angles.  Vortices  further  out 
generate  significant  levels  of  negative  v! ,  but  generate  low  levels  of  positive  v  (§3.1). 
These  vortices  (with  lower  turning  angles)  will  therefore  contribute  more  to  the 
intensity  angle  than  to  the  shear  stress  angle.  Thus,  the  intensity  angle  is  lower 
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because  it  has  more  contributions  from  the  slightly  turned  vortices  further  from  the 
wall  than  the  shear  stress  angle. 

Finally,  a  change  which  is  clearly  seen  in  the  instantaneous  flow  fields  is  the 
kinking  and  later  break-up  of  the  wall  layer  streaks  (Figure  4.36).  It  is  not  known  to 
what  degree  this  phenomenon  affects  the  Reynolds  stresses.  The  break-up  mechanism 
is  primarily  in  the  horizontal  plane.  Vortices  which  in  2D  are  aligned  with  each  other 
in  the  streamwise  direction  next  to  a  streak  turn  in  3D.  They  will  then  no  longer  be 
aligned  behind  each  other  and  the  streak  breaks  up.  Since  intense  high-speed  regions 
on  the  high-speed  streaks  are  well  correlated  with  near-wall  vortices,  these  regions 
turn  swiftly  and  significantly  (Figure  4.36). 

In  the  late  period,  turbulent  intensities  near  the  wall  begin  to  recover.  The  early 
part  of  the  apparent  recovery  is  not  due  to  diminishing  of  3D  effects.  It  will  be 
shown  below  that  the  3D  effects  are  still  present;  however,  they  are  overcome  by 
the  increased  mass  flow.  Recovery  from  3D  effects  should  occur  when  the  vortices 
begin  to  align  with  the  mean  velocity  gradient  direction.  In  order  to  find  when  this 
occurs,  it  is  necessary  to  estimate  the  angle  of  these  vortices  as  a  function  of  time  and 
distance  from  the  wall.  One  way  of  doing  this  is  to  compute  the  direction  of  strongest 
correlation  for  quantities  which  best  capture  vortices.  The  expression  for  this  angle  in 
the  horizontal  plane,  is  given  in  eqn.  (2.7.2).  In  Figure  4.37,  the  angles  of  strongest 
auto-correlation  for  p' ,  V2p'  and  Jx  cos(^)-f  Jz  sin(0#)  are  shown  at  £=0.6.  All  three 
angles  are  nearly  the  same  (except  in  the  viscous  sublayer).  The  mean  velocity  in  the 
direction  normal  to  B %  increases  continuously  in  magnitude  during  the  time-period 
studied  (Figure  4.38).  Thus,  the  three-dimensional  effects  in  the  late  period  have  not 
diminished  and  recovery  is  caused  by  the  increased  “streamwise”  mass  flow. 

4.3  Mean  velocities 

The  mean  velocity  component  in  the  z-direction,  U,  is  virtually  unchanged  during 
the  initial  period  ( t  <  0.3).  Later,  U  is  reduced  near  the  wall  and  since  the  mass  flow 
in  the  z-direction  is  fixed,  it  increases  further  out.  The  percentage  reduction  is  largest 
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very  close  to  the  wall,  where  undergoes  its  maximum  drop  of  about  8%  (Figure 
4.39).  This  is  the  only  observed  quantity  which  was  still  dropping  at  £=1.05,  when 
the  accumulation  of  statistics  was  ended.  During  the  initial  period,  the  weakening 
of  the  events  with  highest  and  lowest  u  cancelled  their  respective  contributions  to 
U,  In  the  later  period,  the  turning  reduces  the  magnitude  of  positive  u'  more  than 
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the  magnitude  of  negative  u'  because  the  vortices  closest  to  the  wall  (which  turn) 
generate  the  most  intense  positive  u1  while  vortices  further  out  (which  don’t  turn) 
generate  as  intense  negative  u'  as  vortices  close  to  the  wall.  The  net  effect  is  a  reduced 
U  near  the  wall  in  the  later  period. 

The  spanwise  mean  velocity  component,  W,  is  identical  to  the  laminar  profile 
(§2.5)  during  the  initial  period.  When  the  structures  turn,  they  make  the  Reynolds 
shear  stress  vector  turn  as  well  and  the  W  profile  is  no  longer  laminar  (Figure  4.40). 
However,  W  does  not  develop  a  logarithmic  region  during  the  time-period  studied. 

4.4  Vorticity  intensities 

In  2D  channel  flows,  the  x-vorticity  intensity,  wj.2,  is  mostly  due  to  quasi- 
streamwise  vortices  and  the  secondary  vorticity  they  create  below  them  (Kim  et  al. 
1987).  Because  the  sign  of  the  secondary  vorticity  is  opposite  to  that  of  the  primary 
vortex,  there  is  a  local  minimum  in  the  viscous  sublayer  (Figure  4.41a).  At  <=0.3,  u42 
is  reduced  all  the  way  out  to  y%D  ~  80  with  the  largest  percentage  reduction  of  9%  at 
the  wall.  The  reason  for  the  reduction  near  the  wall,  where  %  ?//2/j/2,  has  already 
been  discussed  in  Section  4.1.  The  reduction  outside  the  spanwise  boundary  layer 
is  believed  to  be  due  to  reduced  gradients  caused  by  the  weakened  ejections  from 
the  near-wall  region.  After  <=0.3,  the  shear  layers  which  contain  strong  spanwise 
vorticity  begin  to  turn.  Since  is  significantly  larger  than  uJ}  in  near- wall  2D  flows, 
the  turning  tends  to  increase  u'£  in  this  region.  Outside  y£D  —  50,  the  two  vorticity 
intensities  are  of  about  equal  magnitude,  and  wj?  is  reduced  due  to  weakened  ejections 
from  the  near-wall  region.  The  weighted  PDF  for  u)'x  at  y+D  =  20  shows  an  increase 
in  negative  values  and  a  decrease  in  positive  values  at  <= 0.6  (Figure  4.42a).  This  is 
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consistent  with  the  turning  of  the  shear-layers  which  generate  more  <  0.  Note 
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that  the  turning  transforms  some  of  ^-  <  0  of  the  shear  layers  to  <  0. 

Figure  4.416  shows  that  the  vertical  vorticity  intensity,  drops  across  the 
boundary  layer.  Note  that  u>y  is  invariant  with  coordinate  rotation.  The  maximum 
percentage  reduction  is  20%  at  <— 0.T5  and  t/+D= 38;  recovery  starts  at  <=0.5  near 
the  wall.  The  weighed  PDFs  (Figure  4.426)  show  that  both  positive  and  negative 
values  of  significant  magnitude  are  reduced.  These  values  are  usually  found  between 
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near-wall  streaks,  and  next  to  shear  layers  further  out  (high  |^-|  in  both  cases).  The 
reduction  is  therefore  consistent  with  the  reduced  values  of  it' .  Recovery  starts  earlier 
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than  for  the  streamwise  velocity  intensity  and  the  turbulent  kinetic  energy,  which 
indicates  that  the  normal  vorticity  is  more  directly  affected  by  the  break-up  of  the 
streaks  (§4.2.3)  than  the  streamwise  velocity  fluctuations. 

The  spanwise  vorticity  intensity,  (Figure  4.41c)  is  proportional  to  u ®  very 
close  to  the  wall  and  decreases  due  to  the  weakening  and  turning  of  the  streaks. 
During  the  initial  period,  an  increase  occurs  between  y^D  =  7  and  25  because  of  the 
increased  dissipation  associated  with  ejections  generated  by  vortices  with  a  negative 
sign  of  rotation  (§4.1.2).  The  PDF  for  i<4  at  f=0.3  (Figure  4.42c)  shows  an  increase 
in  contribution  from  <JZ  >  0  which  is  consistent  with  more  intense  <  0.  Later, 
this  increase  develops  into  a  plateau  due  to  the  turning  of  the  vortices  (see  profile  of 
<jj'£  Figure  4.41a). 

4.5  Two-point  correlations  and  energy  spectra 

Two- point  velocity  correlations 

Rij(Ax,Az,y,y,t)  =  u'{(x,  z,y,t)u'j(x  +  Ax,z  4-  A  z,y',t)  (4.5.1) 

contain  information  about  the  size  and  relative  locations  of  structures.  In  a  2D  near¬ 
wall  layer,  R\\  is  used  to  study  the  properties  of  the  streaks  while  R22  and  R33  are 
used  to  extract  information  on  the  average  size  and  intensity  of  the  quasi-streamwise 
vortices.  In  3D,  correlation  contour  lines  in  vertical  planes  are  skewed  (shown  for  R\\ 
in  Figure  4.43),  while  in  the  horizontal  plane  they  reflect  the  weakening  and  turning 
of  the  turbulence  structures  (see  below).  Below,  the  correlations  of  the  velocities 
in  the  rotated  coordinates  in  horizontal  planes  are  discussed.  The  rotation  angle 
corresponds  to  the  direction  of  maximum  streamwise  velocity  intensity  (eqn.  2.1  d). 
Thus,  /?,•,*  is  defined  as  in  eqn.  (4.6.1)  with  replaced  by  (§2.1). 

Figures  4.44,  4.45,  4.46  show  the  evolution  of  velocity  two-point  correlations  at 
vto  —  10-  At  t=0.3,  W  =  0.99WOO  at  =  20  and  the  structures  have  turned  very 
little.  The  changes  in  the  correlation  length  and  magnitude  are  therefore  mainly  due 
to  the  weakening  of  the  velocity  fluctuations.  Later,  the  most  pronounced  turning- 
effect  is  the  formation  of  a  kink  in  the  contour  lines  for  R\u  and  a  significant  difference 
in  the  directions  of  maximum  correlation  for  short  separations  between  the  three 
auto-correlations.  The  angle  of  maximum  correlation  for  /?n*  remains  fairly  close 
to  the  Reynolds  shear  stress  angle  for  a  wide  range  of  distances  from  the  wall  and 
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times,  while  the  angles  for  R22  and  #33*  are  larger.  Figure  4.47  shows  the  angles  of 
maximum  correlation  obtained  from  eqn.  (2.7.2)  as  a  function  of  the  distance  from 
the  wall. 

Differences  in  correlation  angles  are  believed  to  be  due  to  the  variation  of  the 
turning  angles  of  the  vortices  with  distance  from  the  wall.  A  review  of  the  discussion 
of  Figure  3.1  regarding  the  relation  between  the  distance  of  a  vortex  from  the  wall 
and  the  velocity  fluctuations  it  generates  in  2D  flows  is  therefore  useful.  The  vortices 
which  are  closest  to  the  wall  turn  the  most  and  are  the  primary  contributors  to  all  the 
two-point  correlations  near  the  wall.  Vortices  further  out  generate  fairly  low  velocity 
fluctuations  near  the  wall  except  for  negative  u*  and,  therefore,  contribute  with  their 
low  turning  angles,  more  to  R{\*  than  to  R22  and  R33*.  The  angle  of  maximum 
correlation  for  is  therefore  the  lowest  of  the  three  near  the  wall. 

The  kink  in  is  believed  to  be  due  to  the  turning  of  the  near  wall  vortices. 
For  large  separations  Ax  and  A 2,  u'(x,  z )  and  u'(x  -f  Ax,  2  +  A2)  are  generated  by 
different  vortices.  Therefore,  the  correlation  angle  tends  to  be  lower  than  the  angle 
for  short  separations  where  the  two  points  would  be  located  on  a  single  vortex  (see 
sketch  Figure  4.48).  R22  and  #33*  show  a  weaker  kink.  This  is  because  the  vortices 
in  Figure  4.48  contribute  to  R22  and  #33*  over  long  distances  only  if  they  both  are 
located  near  the  wall.  On  the  other  hand,  vortices  both  near  and  away  from  the  wall 
contribute  to  intense  u*  <  0  near  the  wall. 

The  two-dimensional  energy  spectra  for  the  u'^  velocities,  Eu*,  E22  and  £33* 
are  shown  in  Figures  4.49.  E\\+  in  the  3D  flow  shows  that  the  large  scales  (low 
frequencies)  have  rotated  less  than  the  the  small  scales.  This  is  another  manifestation 
of  the  kink  in  Ru+-  Ru*  at  large  separations  is  mostly  made  up  of  large  scale 
structures  which  have  turned  by  a  smaller  angle.  The  lack  of  pronounced  kinks  in 
R22  and  /?33*  is  related  to  a  more  uniform  angle  of  spectral  contours  for  E22  and 
#33*- 

4.6  Comparison  with  previously  suggested  mechanisms 

The  mechanisms  for  the  reduction  of  Reynolds  stresses  in  3D  flows  which  have 
been  suggested  in  the  literature  were  described  in  §1.4.4.  These  mechanisms  have 
been  tested  in  the  channel  flow  of  the  present  study.  The  results  of  these  tests  are 
described  below. 
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Lohmann  (1976)  suggested  that  in  a  flow  along  a  cylinder,  part  of  which  is 
spinning,  the  vortices  would  be  torn  apart  when  they  encounter  the  cross-flow.  Such 
a  break-up  was  not  observed  in  the  current  work  (compare  x-vorticity  contours  in 
Figures  4.5  and  4.28),  but  the  vortices  in  Lohmann’s  flow  experienced  a  stronger  and 
more  sudden  cross-flow  which  may  have  caused  additional  effects. 

Vortices  outside  the  inner  layer  tend  to  have  arch  shaped  heads  which  induce 
ejections  (Robinson,  1991).  Bradshaw  and  Pontikos  (1985)  suggested  that  in  3D 
boundary  layers  the  heads  are  tilted  sideways  and  become  less  efficient  in  producing 
Reynolds  shear  stress.  This  hypothesis  can  not  be  tested  in  the  current  flow  due  to 
the  lack  of  inviscid  skewing  mechanism  (see  §1.4.2). 

Anderson  and  Eaton  (1989)  suggested  that  the  mean  streamwise  vorticity,  fl*, 
generated  in  3D,  leads  to  the  cancellation  of  the  natural  boundary  layer  vortices  of 
opposite  sign,  and  also  rolls  up  into  streamwise  vortices  of  the  same  sign.  Thus,  a 
preference  for  vortices  of  one  sign  is  developed,  resulting  in  fewer  strong  ejections  in 
between  vortices  of  opposite  signs.  There  are  two  ways  in  which  ttx  might  reduce 
the  number  of  vortices  of  opposite  sign:  either  the  generation  of  these  vortices  is 
inhibited  and/or  the  existing  vortices  are  destroyed.  The  mechanics  of  the  formation 
of  streamwise  vortices  were  discussed  in  §3.1  for  2D  boundary  layers.  It  was  found 
that  sheets  of  u>x  were  formed  close  to  the  wall  which  rolled  up  into  vortices.  The 
vorticity  of  the  sheet  was  of  opposite  sign  to  the  primary  vortex  that  generated  the 
sheet.  This  process  was  not  significantly  affected  when  a  spanwise  pressure  gradient 
was  imposed.  Using  the  minimal  channel  (§2.4)  as  the  basic  test  set-up,  a  2D  flow' 
with  a  strong  streamwise  vortex  of  negative  sign  was  selected.  Then  two  simulations 
with  the  imposed  pressure  gradient  in  the  positive  and  negative  2-directions  were  run. 
Comparing  the  3  simulations  (2D  and  two  3D’s)  did  not  reveal  a  marked  difference 
in  the  process  of  formation  and  roll-up  of  a  positive  vortex  sheet  into  a  streamwise 
vortex  (Figure  4.50).  Note  that  in  Figure  4.50c  the  vorticity  of  the  vortex  sheet  is 
stronger  than  that  in  2D  because  it  is  superimposed  on  the  mean  streamwise  vorticity. 
However,  the  Figures  showing  the  roll-up  process  (b,  d,  f)  further  downstream  are  very 
similar. 

In  addition,  conditional  averaging  was  used  to  study  whether  there  were  fewer 
vortices  of  one  sign  in  the  3D  channel.  Several  conditions  were  attempted  (various 
ranges  of  w'x  and  v ')  and,  depending  on  the  condition  used,  the  results  sometimes 
indicated  a  weak  preference  for  negative  vortices,  but  never  a  preference  for  positive 
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vortices  (mean  streamwise  vorticity  is  positive  in  the  present  flow).  The  mechanism 
proposed  by  Anderson  and  Eaton  (1989)  therefore  does  not  to  occur  in  this  flow. 

Simpson  and  Devenport  (1990)  argued  that  low-speed  fluid  which  is  pushed 
sideways  by  high-speed  fluid  near  the  wall  need  not  be  ejected  away  from  the  wall  in 
3D  flows  due  to  the  extra  “spatial  degree  of  freedom”.  The  reduction  in  strength  of 
the  ejections  is  consistent  with  the  findings  of  the  present  study  (Figure  4.46  and  c). 
However,  their  reasoning  in  terms  of  a  third  dimension  of  spatial  freedom  is  unclear 
since  the  “three-dimensionality”  only  referred  to  in  the  mean  flow  and  both  2D  and 
3D  flows  are  three-dimensional. 

Eaton  (1992)  suggested  that  the  near-wall  streamwise  vortices  respond  to  3D 
effects  in  the  same  manner  as  the  large  imbedded  vortices  investigated  by  Shizawa 
and  Eaton  (1990)  would.  They  found  that  vortices  with  sign  of  rotation  opposite  to 
♦'■■at  of  the  mean  streamwise  vorticity  interacted  with  the  cross-flow  in  a  very  different 
manner  than  those  of  opposite  sign.  Much  weaker  negative  streamwise  fluctuations 
were  produced  (weaker  ejections)  if  the  vortex- generated  motion  near  the  wall  was  in 
the  same  direction  as  the  mean  cross-flow  while  vortices  of  opposite  sign  generated 
stronger  ejections  of  low-speed  fluid.  The  imbedded  vortices  were  well  outside  the 
location  of  maximum  spanwise  cross-flow  such  that  the  cross-flow  below  them  was 
the  same  in  both  fixed  coordinates  and  coordinates  moving  with  the  vortices  (see 
sketch  in  Figure  4.51).  This  is  not  the  case  for  the  naturally  occurring  streamwise 
vortices  closest  to  the  wall  in  the  3D  channel  flow.  That  is,  in  the  3D  channel  flow 
in  coordinates  moving  with  the  vortices,  the  ejections  are  most  reduced  for  vortices 
with  induced  motion  near  the  wail  opposite  to  the  mean  spanwise  flow  (Figure  4.46). 
It  therefore  seems  that  the  imbedded  vortices  behave  differently  from  the  near-wall 
vortices  in  TBLs.  We  note  that  the  embedded  vortices  of  Shizawa  and  Eaton  were 
much  stronger  than  naturally  occurring  vortices  in  TBLs.  fn  order  to  study  the  effect 
of  the  cross-flow  on  the  most  intense  vortices  in  the  channel  flow,  the  average  flow' 
fields  around  a  strong  vortex  was  computed  using  conditional  averaging.  At  i— 0.0, 
0.3  and  0.6  the  ensemble  averaged  flow  field  around  points  with  |-^~ -|  >  100  at  the 
wall  was  computed  (Figure  4.52).  It  can  be  seen  that  the  conditionally  averaged  flow 
is  consistent  with  the  experimental  result  for  the  imbedded  vortices.  That  is,  strong 
vortices  with  induced  motion  in  the  direction  of  the  mean  spanwise  flow  generate  the 
weakest  ejections.  A  plausible  explanation  is  that  negative  vortices  of  medium  or  low 
strength  are  unable  to  oppose  the  mean  cross-flow  and  lift  the  fluid  closest  to  the  wall 
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as  described  in  §4.1.2  (Figure  iA.b).  However,  the  strongest  ■ ’or t ices  are  able  to  lift 
fluid  from  very  close  to  the  wall  as  well  as  the  fluid  being  converted  toward  them. 

Finally,  Littell  and  Eaton  (1991)  measured  turbulence  statistics  in  an  equilibrium 
3D  turoulent  boundary  layer  on  a  spinning  disk  (see  §1.4.3).  Conditionally  averaged 
two-point  correlations  offered  information  consistent  with  the  present  findings  on 
the  effect  of  the  cross-flow  on  the  streamwise  vortices.  Their  results  indicated  that 
vortices  with  the  same  sign  as  the  local  spanwise  mean  velocity  gradient  generated 
weaker  sweep-motion  while  vortices  of  opposite  sign  generated  weaker  ejections. 
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Figure  4.1.  Reynolds  shear  stress  profiles  in  (a)  initial  coordinates  and  (6) 
coordinates  aligned  with  the  Reynolds  shear  stress  direction  (see  eqn.  2.16).  Starting 
from  i=0.  ( - )  with  increments,  A t  =  0.15. 
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Figure  4.4.  Schematics  showing  the  underlying  mechanisms  I-IV  which  lead  to  the 
reduction  of  the  Reynolds  stresses.  The  observer  is  moving  in  the  spanwise  direction 
with  the  center  of  the  streamwise  vortex. 
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Figure  4.5.  Top  view  of  streaks  (|^|-|  >  50  at  the  wall,  shaded  grey)  and  vortices 

(|u4l  >  60ur2o/^)  for  the  companion  flows  at  t= 0.3.  (a)  2D  flow  and  ( b )  3D  flow. 
Tic  marks  are  0.5  apart. 
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Figure  4.6.  Percentage  reduction  of  streamwise  velocity  intensity  in  (a)  initial 
coordinates  and  (6)  coordinates  aligned  with  the  direction  of  maximum  velocity 
intensity.  Starting  with  <=0.3  with  increments,  At  —  0.15. 
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Figure  4.7.  Percentage  reduction  of  the  normal  velocity  intensity.  Starting  from 
<=0.3  with  increments,  A t  —  0.15. 
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Figure  4.8.  Percentage  change  in  spanwise  velocity  intensity  in  initial  coordinates. 
Starting  with  <=0.3  with  increments,  A t  =  0.15. 


Au*i>' / u*t>, 2Dsl00  IOQAuV/uVjd 


0  20  40  60  80  100  12 0  140 

2/ 2D 


Figure  4.9.  Percentage  reduction  of  streamwise  Reynolds  shear  stress  in  (a)  initial 
coordinates  and  (6)  coordinates  aligned  with  the  Reynolds  shear  stress  direction  (eqn. 
2.16).  Starting  with  t=0.3  with  increments,  A<  =  0.15. 
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Figure  4.10.  Intensity-weighted  PDF  for  (a)  u  at  y?D  =  5,  (6)  v'  at  y+D  =  10  and  (c) 

at  the  wall.  { - ),2=0.;  ( - ),  2=0.3;  ( .  ),  2=0.6.  ?/*  and  w*  are  rotated 

with  the  direction  of  maximum  intensity  (eqn.  2.1  d)  for  2=0.6  and  is  unrotated  for 


Figure  4.19.  A  spanwise  trace  of  u 1  in  2D  (— )  and  3D  (-  -  -)  companion  flows  (a) 
unshifted  and  (b)  shifted  in  the  z-direction. 
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Figure  4.22.  Sample  vortex  from  the  companion  flows  with  positive  sign  of  rotation 

for  (a)  2D  flow  and  (b)  3D  flow  at  <  =0.3.  Contour  lines  below  y?D  —  10  :  (  -  ) 

u'  >  5  with  increments  A it  —  O.o  and  ( - )  ur  <  —  1  with  A u'  —  —0.5.  The  tic 

marks  are  separated  by  0  05. 
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Figure  4.23.  Sample  vortex  from  the  companion  flows  with  negative  sign  of  rotation 
for  (a)  2D  flow  and  ( b )  3D  flow  at  t= 0.3.  Contour  lines  below  y \D  =  10  :  (  - ) 

u'  >  4  with  increments  At/  =  0.5  and  ( - )  v!  <  -1.5  with  At/  =  -0.5.  The  tic 

marks  are  separated  by  0.05. 
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Figure  4.25.  Sample  vortex  with  a  positive  sign  of  rotation  from  (a)  the  2D  and  ( b ) 

r\  / 

the  3D  companion  flows  at  at  t=0.3.  Solid  lines  show  contours  of  <  —100  with 

increment  10.  Dotted  lines  show  u  <  7  with  increment  1.  Tic  marks  are  separated 
yy  0.05. 


FfGURE  4.26.  Schematics  showing  the  events  in  which  u'v'  is  most  reduced  at  y%D  =  25. 
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Figure  4.28.  Top  view  of  streaks  (|^|  >  50  at  the  wall,  shaded  grey)  and  vortices 
(\u)'x\  >  60)  for  (a)  <=0.45  and  (6)  <=0.6.  Tic  marks  are  separated  by  0.5. 
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Figure  4.29.  Angle  of  maximum  Reynolds  shear  stress  ( - )  and  angle  of  maximum 

velocity  intensity  ( - ).  Starting  from  <=0.15  with  increments.  At  =  0.15. 
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Figure  4.30.  Change  immagnitudes  of  the  contributions  from  each  quadrant  for  uW 
in  coordinates  aligned  with  the  Reynolds  shear  stress  at  t=0.6.  u  >  0  and  v'  >  0 
( - ),  u'  <  0  and  v'  >  0  ( - ),  v!  <  0  and  v'  <  0  ( - ),  u'  >  0  and  t/  <  0 
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Figure  4.31.  Conditionally  averaged  field  of  u*  at  yfD  =  5  around  20  <  <  50 

in  a,b,c  and  -50  <  <  -20  in  d,e,f  at  the  wall.  *=(0.,  0.3,  0.6)  in  ( a,b,c )  and 

( d,e,f)  respectively.  The  positive  contours  ( - )  start  from  0.5  with  increment  0.1 

and  negative  contours  ( - )  from  -0.5  with  increment  -0.1.  u*,  te*  are  unrotated 

for  f=0.0,  0.3  and  rotated  witn  the  intensity  angle  (eqn.  2.1  d)  far  t=0.6.  Tic  marks 
are  separated  by  0.1. 
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Figure  4.32.  Top  view  of  the  average  field  of  negative  u*  (  rotated  with  the  intensity 
angle,  eqn.  2.1  d)  at  y^D  =  5  around  -0.25  <  v  <  -0.2  at  y*D  —  10.  Starting  with 
tx*  =  —0.2,  the  increment  is  /W*  =  —0.1.  (a)  t=0.0,  (6)  t—  0.3.  Tic  marks  are  0.1 
apart,  and  (c)  *=0.6, 
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Figure  4.33.  Change  in  magnitude  of  the  contribution  from  each  octant  to  u* v' 
at  <=0.6  in  coordinates  aligned  with  the  Reynolds  shear  stress,  u*  >  0  and  v1  >  0 

( - ),  u*  <  0  and  v'  >  0  ( - ),  u*  <  0  and  v'  <  0  ( - ),  u',  >  0  and  v'  <  0 

( .  ).  Symbols  (x)  represents  tn*  <  0  and  no  symbols  represents  >  0. 
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Figure  4.34.  Regions  of  u'  >  0  (grey)  extend  toward  the  ejection-region  in  3D.  This 
effect  is  weak  during  the  initial  period  as  seen  by  comparing  the  2D  (a)  and  the  3D  ( b ) 
companion  flows  at  t=0.3.  At  t= 0.6  (c)  the  effect  is  much  stronger.  In  addition,  the 
particle  trajectories  will  be  altered  by  the  cross-flow.  This  effect  is  seen  in  (a)  where 

some  “streamlines”  for  both  the  2D  ( - )  and  the  3D  (• - )  companion  flows  are 

shown  traced  from  common  starting  points.  This  effect  will  increase  in  importance 
with  increasing  spanwise  boundary  layer  thickness. 
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the  average  vortex-angle  at  y?D  —  20  and  move  with  the  mean  flow  at  this  location 
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Figure  4  39.  Percentage  reduction  in  the  x-component  of  the  mean  velocity.  1'. 
starting  from  t=0.15  with  increments,  At  =  0.15. 
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Figure  4.40.  Evolution  of  laminar  ( - )  and  turbulent  (  -  )  spanwise  mean 

velocity  profiles,  W ,  starting  from  <=0.15  with  increments,  A<  =  0.15. 


Figure  4.43.  R\\  at  (a)  t—0.  and  ( b )  t— 0.6  around  y%D  =  10.  Positive  contour 

levels  are  incremented  by  0.5,  starting  from  1  (  -  )  and  negative  contours  are 

incremented  by  —0.5,  starting  from  —0.5  ( - ).  Tic  marks  are  0.1  apart. 


Figure  4.44.  /?u*  (based  on  the  velocity  components  in  the  direction  of  maximum 

intensity)  at  y^n  =  10  for  (a)  1=0.0,  (6)  1=0.3  and  (c)  t=0.6.  Positive  contour  levels 

are  incremented  by  1,  starting  from  1  ( - )  and  negative  contours  are  incremented 

by  —0.25,  starting  from  —0.5  ( - ).  Tic  marks  are  0.1  apart. 
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Figure  4.45.  R22  a*<  vto  =  10  for  (a)  <=0.0,  (6)  <=0.3  and  (c)  <=0.6.  Positive  contour 

levels  are  incremented  by  0.01,  starting  from  0.01  ( - )  and  negative  contours  are 

incremented  by  —0.01,  starting  from  —0.01  ( - ).  Tic  marks  are  0.1  apart. 


Figure  4.46.  Rzz*  (based  on  the  velocity  components  in  the  direction  of  maximum 
intensity)  at  y^D  =  10  for  (a)  <=0.0,  ( b )  t— 0.3  and  (c)  <=0.6.  Positive  contour 

levels  are  incremented  by  0.1,  starting  from  0.2  ( - )  and  negative  contours  are 

incremented  by  —0.005,  starting  from  —0.025  ( - ).  Tic  marks  are  0.1  apart. 


Figure  4.48.  View  of  two  vortices  from  below  (grey)  in  3D  flow.  Angle  for  short 
streamwise  separation  (0i)  is  larger  than  the  angle  for  large  separations  (#2)- 
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Figure  4.50.  Formation  of  a  vortex  sheet  with  ux  >  0  in  flow  with  (lx  —  0  (a,  6), 
nx  >  0  (c,  d),  Ox  <  0  (e,  J).  Figures  to  the  left  show  plane  across  aft  of  the  sheet 
and  figures  to  the  right  show  a  plane  across  the  front  part  of  the  sheet.  Contour  lines 

represents  u>x  with  increment  0.25  starting  with  0.5  ( - )  and  with  increment  —0.25 

starting  with  —0.5  ( - ).  Grey  regions  show  areas  with  u>x  <  0. 


109 


(a) 


(b) 


Figure  4.51.  Schematics  showing  the  spanwise  mean  velocity  profile  in  coordinates 
aligned  and  moving  with  a  vortex  which  is  (a)  above  and  (b)  near  the  location  of 
maximum  cross-flow.  Even  though  the  profile  is  the  same  in  coordinates  fixed  to  the 
wall,  the  spanwise  mean  velocity  below  the  two  vortices  are  of  opposite  sign. 
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Figure  4.52.  Top  view  of  average  field  of  u*  (aligned  with  maximum  streamwise 
intensity)  at  y%D  =  5  around  >  100 uT/S  (a,  b,  c)  and  <  100ur/6  (d,  e,  f).  t 

=  0.,  0.3  and  0.6  (top  to  bottom).  The  increment  is  0.5  starting  with  1.5  ( - )  and 

—0.25  starting  with  —0.5  ( - ). 
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CHAPTER  5 


CONCLUSIONS 

T.  e  effect  of  an  impulsively  started  transverse  pressure  gradient  on  turbulence 
structures  in  a  planar  channel  flow  was  studied  using  direc'  numerical  simulation. 
In  agreement  with  experiments,  the  simulation  showed  a  reduction  in  the  Reynolds 
shear  stress.  It  also  showed  a  drop  in  the  turbulent  kinetic  energy.  To  understand  the 
underlying  mechanisms  for  these  seemingly  counter-intuitive  observations  we  focused 
our  study  on  the  effects  of  three-dimensionality  on  the  wall  layer  structures.  To 
this  end  we  performed  a  detailed  study  of  the  dynamics  of  the  energy  producing 
structures  in  both  2D  and  3D  channel  flows.  This  study  included  detailed  examination 
of  the  instantaneous  flow  fields,  probability  density  functions,  two-point  statistics 
and  conditional  averaged  fields.  The  mechanisms  described  are  consistent  with  all 
statistical  data  as  well  as  observations  of  instantaneous  fields. 

In  the  2D  channel  flow  our  study  of  the  near  wall  structures  revealed  a  mechanism 
for  the  generation  of  streamwise  vortices.  The  cycle  of  events  begins  as  follows.  An 
inclined  sheet  of  streamwise  vorticity  is  generated  below  an  existing  inclined  stream- 
wise  vortex  through  tilting  of  normal  vorticity.  The  sheet  is  unstable  and  rolls  up  into 
an  inclined  streamwise  vortex  which  will  grow  due  to  stretching.  The  key  element 
for  the  stretching  is  the  inclination  of  the  vortex  and  the  transport  of  fluid  by  the 
vortex  in  the  wall-normal  direction.  The  extent  of  stieamwise  vortices  in  the  channel 
were  rarely  more  than  Ax+  =  350.  Near-wall  streaks  are  generated  and  maintained 
by  a  string  of  these  vortices.  Low-speed  streaks  are  longer  than  high-speed  streaks 
because  the  magnitude  of  negative  u'  generated  by  a  vortex  is  fairly  independent  of 
the  distance  of  the  vortex  from  the  wall  and  her.ee  von'ces  at  various  locations  con¬ 
tribute  to  a  low  speed  streak.  On  the  other  hand  an  increase  in  the  height  of  a  vorcex 
above  the  wall  tends  to  decrease  the  magnitude  of  positive  u'  it  generates  near  the 
wall. 

In  the  3D  channel  the  reduction  in  total  Reynolds  shear  stress  and  turbulent  ki¬ 
netic  energy  were  shown  to  be  due  to  changes  in  the  particle  trajectories  associated 
with  streamwise  vortices.  Most  affected  are  trajectories  on  the  sweep-side  of  vortices 
which  move  fluid  below  them  in  the  same  direction  as  the  spanwise  mean  velocity  in 
coordinates  moving  with  the  vortex,  and  on  the  ejection-side  of  vortices  of  opposite 
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sign.  These  sweeps  are  weakened  because  the  fluid  moving  toward  the  wall  is  con- 
vected  toward  the  vortex-center  where  the  vertical  vortex- induced  motion  is  low.  The 
fluid  will  therefore  not  reach  as  close  to  the  wall  in  3D  and  hence  its  contribution 
to  the  Reynolds  stresses  is  reduced.  The  Reynolds  stresses  generated  by  ejections 
due  to  vortices  of  opposite  sign  are  reduced  because  the  lifted  fluid  originates  further 
away  from  the  wall  in  3D.  In  addition  these  ejections  are  reduced  due  to  increased 
viscous  dissipation  when  low-speed  fluid  is  convected  sideways  above  high-speed  fluid. 
Reynolds  stresses  generated  by  other  sweep  and  ejection  motions  are  also  reduced  due 
to  changes  in  particle  trajectories. 

Among  the  other  observed  phenomena  of  interest  were  the  break-up  of  the  wall 
layer  streaks  during  the  later  period  of  three-dimensionality.  This  break-up  was  at¬ 
tributed  to  the  differential  turning  of  the  streamwise  vortices  according  to  their  dis¬ 
tance  from  the  wall  and  with  respect  to  the  streaks. 

In  summary,  our  numerical  study  of  the  various  stages  of  development  of  a  3D 
flow  driven  by  a  spatially  homogeneous  spanwise  pressure  gradient  has  illuminated  the 
reasons  why  the  extra  increase  in  the  total  shear  rate  leads  to  decreases  in  turbulent 
kinetic  energy  and  total  shear  stress.  Future  work  should  include  experimental  veri¬ 
fication  of  this  study.  In  most  cases  experimental  verification  is  feasible  only  through 
novel  conditional  averages  and  two-point  correlations.  Direct  numerical  simulations 
of  a  spatially  developing  3D  turbulent  boundary  layer  is  highly  desirable.  Such  a 
computation  will  include  inviscid  skewing  effects  which  is  an  important  feature  of 
many  of  the  experimentally  studied  3D  flows. 
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